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ABSTRACT 

vn 

We present here a new spectral energy distribution (SED) fitting approach that we adopt to select radio-excess sources amongst distant 
^— ^ star-forming galaxies in the GOODS-Herschel (North) field and to reveal the presence of hidden, highly obscured AGN. Through 

I extensive SED analysis of 458 galaxies with radio 1.4 GHz and mid-IR 24 /jm detections using some of the deepest Chandra X-ray, 

Spitzer and Herschel infrared, and VLA radio data available to date, we have robustly identified a sample of 5 1 radio-excess AGN 
I (~1300 deg"-) out to redshift z ~ 3. These radio-excess AGN have a significantly lower far-IR/radio ratio (q < 1.68, 3cr) than the 

• . typical relation observed for star-forming galaxies {q x 2.2). We find that «45% of these radio-excess sources have a dominant AGN 

^ component in the mid-IR band, while for the remainders the excess radio emission is the only indicator of AGN activity. The presence 

of an AGN is also confirmed by the detection of a compact radio core in deep VLBI 1 .4 GHz observations for eight of our radio-excess 
rN sources (a 16%; a;66% of the VLBI detected sources in this field), with the excess radio flux measured from our SED analysis agreeing, 

to within a factor of two, with the radio core emission measured by VLBI. We find that the fraction of radio-excess AGN increases 
with X-ray luminosity reaching ~60% at Lx ~ IG^ - 10"*^ erg s"', making these sources an important part of the total AGN population. 
However, almost half (24/51) of these radio-excess AGN are not detected in the deep Chandra X-ray data, suggesting that some of 
these sources might be heavily obscured. Amongst the radio-excess AGN we can distinguish three groups of objects: i) AGN clearly 
identified in infrared (and often in X-rays), a fraction of which are likely to be distant Compton-thick AGN, ii) moderate luminosity 
AGN (Lx g lO""' erg s"') hosted in strong star-forming galaxies, and iii) a small fraction of low accretion-rate AGN hosted in passive 
(i.e. weak or no star-forming) galaxies. We also find that the specific star formation rates (sSFRs) of the radio-excess AGN are on 
average lower that those observed for X-ray selected AGN hosts, indicating that our sources are forming stars more slowly than typical 
AGN hosts, and possibly their star formation is progressively quenching. 

Key words, galaxies: active - quasars: general - X-rays: galaxies - infrared: galaxies - galaxies: star formation 
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1. Introduction 

The discovery of a strong correlation between the properties 
of galaxies and those of the supermassive black holes (SMBH) 
hosted in their centres, such as t he Mbh - Mbuige or Mb h - cr 
tions (FeiTarese & Merritt'2000', 'Gebha rdt et al.|2 000; MagoiTian] 



etaLll998, Marconi & Hunt 2003), has pointed out that SMBHs 



must play an important role in the growth and evolution of galax- 
ies (see Alexander & HickoX|2012 for a general review). In the 
past decades many studies have focussed on understanding the 
relation between nuclear activity (AGN) and host galaxies and 
have revealed a common history, where both star formation and 
black hole accretion were much more common in the past, with a 
peak at redshift z * 2 (e.g. 'Fiore et al.|2003| IHopkins et al.|20071 
p006bi Marconi et al. 2004; Merl oni et aL |'2004 2007Jl. Com- 
pleting the census of AGN activity, especially at redshifts where 
most of the accretion occurred, is therefore essential in order to 
understand the nature of the link between SMBH and galaxies 
and their cosmic co-evolution. This paper aims to expand our 
knowledge of the AGN population, by selecting objects with ra- 
dio emission in excess of that expected from star formation. As 
discussed below, this method selects many AGN that cannot be 
identified using other established techniques, and so moves us 
closer to a complete census of growing SMBHs in the Universe. 

Deep X-ray surveys have proved to be a very powerful tool 
in detecting obscured and unobscured AGN down to faint fluxes 
and to high redshifts (z ~ 5; e.g. iAlexand er et al.||2001[ 



etaLf2 009; Fio re et al.|2003tpasm ger 2008l, |LehmeretaL 



Brusa| 



2012 



Xue et al.|201ii However, it is now evident that even the deep 
est X-ray surveys are not complete (e.g. H asinger||2008[ |Tozzi| 
|et al.|200^ and miss a significant part of the AGN population, in 
particular the most obscured, Compton-thick (CT) AGNs, where 
the X-ray emission below 10 keV is strongly suppressed by large 
column density gas (A^h > 10^** cm"^). A large population of 
heavily obscured AGN is indeed predicted by synthesis models 
of the X-ray background (XRB; Ballantyne et al. 201 1[ Comastri] 
let al.|1995i[Gini et al.|2007ir2001t|TreTster et ai|2009l ) in order 
to reproduce the high energy peak of the observed X-ray back- 
ground emission (E ^ 30 keV), which has not yet been directly 
resolved by current X-ray surveys. 

Since large amounts of gas and dust are responsible for the 
suppression of the radiation in the UV, optical and soft X-ray 
bands, perhaps the most obvious waveband to search for these 
heavily obscured objects is where the dust emission peaks, i.e. 
the infrared (IR) band. In fact, the dust surrounding the SMBH 
is heated by the nuclear radiation, reaching temperatures T ~ 
200 - 1000 K, and re-emits the radiation predominantly in the 
mid-infrared (MIR; X * 5-40 jum) band, peaking at ~ 20-30 //m 
(Muflane y et al.||201l} [Netzer et aL]|2007| l. Moreover, at these 
wavelengths the effects of extinction are small, making it eas- 



ier, in theory, to find even the most obscured AGN (e.g. Gandhi 
et al. 1 2009 ; Gouldi ng et al.|2012| l. The downside of using the IR 



band to search for AGN activity is that dust is present not only 
in the circumnuclear region of AGN, but also in the host galaxy, 
in particular in star-forming regions. The dust in these regions is 
heated on average to lower temperatures, and therefore its emis- 
sion peaks at longer wavelengths, than that around the black hole 
(typically at ^ ^100 jim, T j,,,, - 20- 50 K; e.g. [Calzetti et aT 
|2000[[Chary & Elbaz|2001| l. However star formation often dom- 
inates the spectral energy distribution (SED) over the entire IR 
band (e.g. Elvis et al.| 1994[ [Mullaney et"aL][20n] [Netzer et al.] 
^007 i Richards et al. 2006|l and it is often not trivial to separate 
it from the AGN emission. 

Potentially a very powerful approach to obtaining an unbi- 
ased look at obscured and unobscured AGN is through radio 



observations. In fact, at radio frequencies, where the emission 
is mainly due to non-thermal processes, such as synchrotron 
radiation, the effects of extinction are negligible. Historically, 
AGN detected in radio surveys have been divided into two main 
classes: i) radio-loud (RL) AGN, which are the strongest radio 
emitters (typica lly Uad > 10^'^'^^ W Hz [Miller et al.||l990 
Yun et al.|2001] ) and show strong extended radio emission, such 
as kpc-scale relativistic jets and lobes, and ii) radio-quiet (RQ) 
AGN, the weaker radio emitters, whose radio emission is con- 
fined in a small, unresolved region (< 0.1 pc; "core"); the lat- 
ter group constitutes the majority of the population (~90%; e.g. 
[Miller et al.|1990[ [Stocke et al.|1992| l. The separation between 
radio-loud and radio-quiet AGN has typically been set atR - 10, 
where R, the radio-loudness parameter, is defined as the ratio be- 
tween the monochromatic flux density in the radio and optical 
band^/? = 5rad/5opt (e.g. jKellermann et al.| 19891 [Laor ^OOO'- 
However, more recent studies based on deep radio surveys have 
shown that while the radio-loudness parameter spans a very wide 
range of values for AGN, there is no clear evidence of bimodal- 



ity in the population (e.g. 


Brinkmann et al. 2000 Cirasuolo et al. 


[20031 [La Franca et al.|20 


1 0| [White et al.|2000|l. AGN can there- 



fore be identified in the radio band with a wide distribution of 
radio power 

While at bright fluxes the radio population is almost en- 
tirely composed of AGN, at low radio fluxes (sub-mJy regime) 
star-forming galaxies (SFGs) constitute a significant fraction of 
the radio source population and become dominant at fiJy fluxes 
(e.g. [Seymour et al.n2008) . The non-thermal radio continuum 
observed in star-forming galaxies is produced by synchrotron 
radiation from cosmic ray electrons and positrons, accelerated 
by supernova remnants, which mainly occur in young stellar 
populations in star-forming regions (see [Condon[[1992[ for a 
review). The radio emission observed in star-forming galax- 
ies tightly correlates with the emission in the far-infrared (FIR; 
A X 40-120 jum) band, since they both originate from star forma- 
tion processes (e.g. 'Apple ton et al. 2004; Condon 1992; Helou| 
[et al. 1985; Ivison et al. 2010; Yun et al. 2001f7TEs"correlation, 
observed primarily in local star-forming galaxies and starbursts, 
is found to hold out to high redshifts (z ~ 2; e.g. Bourne et al. 
20TT][lbEet al.[2008|[Ivison et al.[2010|[Mao et al.[201 l||Sargent 



et al.|2010b| l 

Joint analyses in the FIR and radio bands, therefore, allow 
us to separate star-forming galaxies from the AGN population. 
Although the weakest radio AGN (RQ) have been found to fol- 
low the same FIR/radio correlation of star-forming galaxies (e.g. 
Moric et al.,2010, Padovani et al.|20TT) , given their wide range 



in radio-loudness, it is possible to identify AGN via their devi- 
ation from the expected FIR/radio relation, the so-called "radio- 
excess" sources (e.g. Donley et al. 2005, Roy & Norris 1997). 
Since the radio emission is not affected (or only lightly affected) 
by extinction, the radio-excess source selection can potentially 
identify AGN that are often missed in optical or even deep X- 
ray surveys (i.e. the most obscured Compton-thick AGN; e.g. 
[Donley et al.l2005). 

In this work, we combine two methods to identify the pres- 
ence of AGN in star-forming galaxies out to high redshift (z ~ 3): 
i) detailed IR SED decomposition, which allows us to measure 
the AGN contribution to the total SED, typically dominated by 
star formation emission; ii) radio-excess signature compared to 
the typical FIR-radio relation observed for SFGs, which is most 



' Most recently, other definitions of the radio-loudness parameter have 
been used, e.g. Rx = log (yLjad/L2-io kcv)> which uses the monochro- 
matic radio luminosity and hard X-ray luminosity (e.g. [B allantyne|2009[ 
La Franca et al.|2010b. 
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likely due to the presence of nuclear activity. We investigate here 
the radio-excess sources in the GOODS-North field, using deep 
infrared, radio and X-ray data, which are some of the deepest 
data available to date. The paper is organised as follows: the data 
and catalogues used in our investigation are presented in Section 
|2] In Section [3] our SED fitting approach is described in details 
as well as the definition of the FIR-radio flux ratio {q) and the 
radio-excess sample selection, together with a comparison with 
other selection criteria used in previous studies. In Section|4] we 
investigate the X-ray, radio and IR properties of our radio-excess 
AGN sample and their SEDs. In Section [5] we discuss the mixed 
population found amongst our radio-excess AGN, attempting to 
constrain the fraction of candidate Compton-thick AGN; we also 
examine the star formation properties of the radio-excess AGN 
hosts through their specific star-formation rate (sSFR) in com- 
parison with those of X-ray selected AGN hosts. In Section [6] 
we summarise our results and give our conclusions. In Appendix 
A, the tests performed to refine our SED fitting approach are ex- 
plained and the best-fit SEDs for the entire radio-excess sample 
are shown in Appendix B. 

Throughout the paper we assume a cosmological mo del with 
Hp = 70 k m s-i Mpc-' 
|etal.|2003| l. 



Q 



M 



0.27 and Qa = 0.73 (Spergel 



noise ratio (S/N) > 3 in the GOODS-N field. However, we note 
that in the outer regions of the GOODS-N field the MIPS data 
is shallower and the uncertainties on the source fluxes are typi- 
cally larger; we therefore limit our catalogues to a smaller area 
(-135 arcmin^) within the GOODS-N field where the MIPS data 
is deeper (1943 sources detected at 24 /zm, ~76%). We require 
at least a detection (S/N > 3) at 24 fim for the sources in our 
sample to be able to constrain the source SEDs in the MIR band 
(see SectirT). 

An area of ~150 arcmin^ of the GOODS-N field has also 
been surveyed at 16 fim using the Infrared Spectrograph (IRS) 
peak-up Imaging (PUI) with pointings of ~10 min each. The ob- 
servations and data reduction are described in detail by T eplitz| 
|et al. (2011 ); the resulting mosaic image is characterised by 0.9" 
pixel scale and has an average Scr depth of ~40 yuJy (Teplitzj 
|et al. 2011). The source catalogue was constructed using Spitzer- 



MIPS 24 jum priors (> Scr sources) and the 16 yum fluxes were 
calculated through PSF-fitting, similarly to the procedure used 
for the 24 jum data; the 16 /im catalogue contains 770 sources 
(Daddi et al. 2012, in prep.). 



2.2. GOODS-Herschel FIR data 



2. Observations and catalogs 

The Great Observatories Origins Deep Survey-North field 
(GOODS-N; [Giavalisco et al.|2004 ^ is one of the deepest multi- 
wavelength surveys currently available and it constitutes an un- 
precedented resource in terms of its broad-band coverage and 
sensitivity. It covers a; 160 arcmin^ centre d on the Hubble Deep 
Field North (HDF-N, 12'' 36'", +62° 14'; [Williams et al.| 1996 1 



and it includes very deep X-ray Chandra data (2 Ms; [Alexander 



et al."2003), optical Hubble Space Telescope (HST; [Giavahsco 



et al. 2004 ) and mid-infrared (MIR) Spitzer observations (PI: M. 
Dickinson); the GOODS-N field has also been the target of sev- 
eral deep optical imaging and spectroscopic campaigns from 8- 
10 m ground-based telescopes. Recently, new deep observations 
of this field in the far-infrared (FIR) band with Herschel (Elbaz] 
let al.|2011| and the radio band with VLA ( [Morrison et al. 2010) 
have usefully increased the potential of the GOODS-N data set. 

2.1. Spitzer MIR data 

The GOODS-N field has been observed at MIR wavelengths 



by Spitzer at 3.6, 4.5, 5.8 and 8.0 yum with IRAC ( [Fazio et alT] 
|2004) , with a mean exposure time per position of x.Ti hours per 
band, and at 24 pm with MIPS (Riek e et al.[[2004| l, as part of 
the GOODS Spitzer Legacy program (PI: M. Dickinson). The 
source catalogue was produced using the SExtractor source de- 
tection routine (Bertin & Arnouts 1996) on a combined 3.6 /jm + 
4.5 fim image, with matched aperture photometry performed in 
the four IRAC bands individually (M. Dickinson et al., in prepa- 
ration). The resulting IRAC catalogue includes 19437 objects 
detected at 3.6 pm with a ~50% completeness limit of 0.5 piy. 

The 24 fim observations consist of a final mosaic image of 
1 .2" pixel scale and a 5cr sensitivity limit of ~30 /iJy. The source 
extraction was performed with a PSF fitting technique using the 
positions of the IRAC 3.6 pm sources detected at > 5cr as priors 
(see Magnelli et al.|201 1 for details). The IRAC 3.6 fim data is 
used to define the source priors because it is ~30 times deeper 
than the 24 pm observations, and therefore all real 24 pm de- 
tected sources should also be detected at 3.6 pm. The resulting 
24 pm catalogue includes 2552 sources detected with signal-to- 



The GOODS-N field has been observed by the Herschel Space 
Observatory as part of the GOODS-Herschel survey (PI: D. El- 
baz), which consists of deep FIR observations of the GOODS- 
North and GOODS-South fields for a total exposure of 361.3 
hours. Imaging of the fufl norther n field (GO ODS-N; 10' x 16') 
was performed using PACS (Poglitsch et al. 2008) at 100 pm and 
160 pm (124 hours of observations) and SPIRE ( [Griffin et al.[ 
2010|l at 250 pm, 350 pm and 500 pm (31.2 hours in total); see 
Elbaz etal.[ ( [20n] ) for details on the PACS and SPIRE observa- 
tions. 

The data reduction was performed following the procedure 
described in Berta et al. ( [201 01) and the resulting images have 
pixel scales of 1.2" and 2!4^at PACS 100 pm and 160 pm, re- 
spectively and 3.6", 5.0" and 7.2" at SPIRE 250 pm, 350 pm and 
500 pm. Since the 350 pm and 500 pm data suffer from strong 
source blending due to the large pixel scales, we did not include 
data at these wavelengths in our SED fitting procedure (Sect 



3.1 



and we only used the PACS 100 pm and 160 pm and SPIRE 250 
pm data. The 350 pm and 500 pm flux densities and/or upper 
limits (S 350 ~ 20.0 mJy and S 500 ~ 30.0 mJy, Scr) from the cat- 
alogue described in Elbaz et al. (2011 ) were only used in some 
of the plots (see e.g. Fig. [B.l[ ) to verify the accuracy of our SED 
fitting at these longer wavelengths. 

The Herschel 100 pm, 160 pm and 250 pm fluxes were calcu- 
lated using PSF fitting at the positions of the Spitzer-MlPS 24 pm 
sources, which are used as priors (Daddi et al. 2012, in prep.). 
In the resulting catalogue 819 sources (~42%) have S/N > 3 in 
at least one of the Herschel bands: 633 sources at 100 pm, 537 
sources at 160 pm and 435 at 250 pm, with 5cr (3cr) sensitivity 
limits of -1.7 (-1.2) mJy, -4.5 (-2.3) mJy and -6.5 (-4.0) mJy, 
respectively. It is important to note, however, that the sensitivity 
of the 250 pm data strongly varies across the field, depending 
on the local source density of the 24 pm priors. We included 
in our sample all of the 1943 24 pm detected sources within the 
GOODS-Heric/ie/ field (restricted to the area with deeper MIPS 
data), with or without a significant Herschel detection; this is to 
avoid biases against faint FIR sources, which are more likely to 
be AGN dominated. 
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2.3. VLA radio data 

Deep, high-resolution radio observations of the GOODS-N field 
were taken at 1.4 GHz using the National Radio Astronomy Ob- 
servatory's (NRAO) Very Large AiTay (VLA) in the A, B, C and 
D configurations (165 hours). The combined radio image reaches 
a rms noise level of ~ 3.9 yuJy beam"' near the centre with a 
beam size of ~ 1.7". These are amongst the deepest radio data 
taken so far. From the VLA image the radio flux density mea- 
surements have been obtained through PSF fitting at each 3.6 //m 
source position (coinciding with the 24 fim positions; see Sect. 
|2.1[ Daddi et al. 2012, in prep.): 1.4 GHz flux measurements 
were obtained for all 1943 sources detected at 24 /urn within the 
GOODS-//er5c/ie/ area, with 489 sources having S/N > 3 and 
Sy > 13 fiiy. The remaining 1454 sources have S/N < 3 and 
they are considered as radio upper limits. Within this field, we 
estimated that ~20 radio detected sources (5cr) are undetected 
at 24 yum down to a flux limit of Sv ~ 21 /iJy; this gives an esti- 
mate of the completeness of our radio and 24 jjm detected sample 
(hereafter VLA/24 fim sample) of ~93%, as compared to a pure 
radio-selected sample. We note that the detection limit adopted 
for our VLA radio catalogue is lower than that used in |Morrison| 
|et al.| ( |2010| l (5cr detection threshold) and therefore the number 
of radio detections found here (489 sources) is much larger than 
that found by |Morrison et al.| ([2010| over the same area (256 
sources). 



the VLA/24 jjm sources, yielding a total of 159 X-ray detected 
sources (i.e. ^33% of the VLA/24 fim sample). We note that 
since the source catalogues in all of the MIR and FIR bands con- 
sidered here, as well as the VLA radio catalogue, are based on 
the 3.6 i-im positions there was no need to cross-match the X- 
ray source positions with any of the other bands. For the X-ray 
undetected sources, 3cr upper limits were derived from aperture- 
corrected photometry in the Chandra images at the 24 yum source 
positions, assuming a power-law model with F = 1.4 (see e.g.. 
Sect. 3.4.1 of ,Alexander et al. (2003,, Bauer et al.|2010[ ). 

The X-ray luminosities (2-10 keV; rest- frame) of the sources 
were extrapolated from the observed 2-8 keV fluxes calculated 
from detailed X-ray spectral analysis (Bauer et al., in prepara- 
tion) and from [Alexander et al.| ( |2003) for the sources in the sup- 
plementary Chandra catalogue. We used the redshifts described 
in Sect. |2.5| and assumed a constant photon index F = 1 .9 (which 
gives a band conversion factor L2-10 keV = 1-08 L2-8 kev)> includ- 
ing appropriate fe-correction. We note that the X-ray luminosi- 
ties have not been corrected for absorption because the column 
density estimates (A^h. available from Bauer et al., in prep.) of- 
ten have large uncertainties; moreover, A^h values were not avail- 
able for the sources in the supplementary catalogue ( Alexander] 
let al.|[2003). Therefore, to avoid adding further uncertainties to 



the X-ray luminosities and to keep consistency in the L2-10 keV 
measurements between the main and the supplementary Chan- 
dra catalogues, we did not apply any absorption corrections. 



2.4. Chandra X-ray data 



2.5. Redshifts 



The Chandra X-ray observations of GOODS-N field cover an 
area of ^ 448 arcmin^ in the 0.5-8.0 keV energy band, with 
an exposure of *2 Ms (Chandra Deep Field North, CDF-N; 
Alexander et al. 2003 | l, reaching a sensitivity (on-axis) of ^ 2.5 x 
10"" ergcm"^ s"' (0.5-2.0 keV) and 1.4 x 10'^^ erg cm ^ s ' 
(2-8 keV). The main source catalogue of the CDF-N includes 
503 X-ray detected sources (Alexande r et al.f2003 



A supple- [2005 



mentary catalogue containing 430 X-ray sources is also available 
in this field. These catalogues where constructed using a source 
detection algorithm with false-positive probability threshold of 
10~^ for the main catalogue, and a mor e relaxed threshold of 
I0~^ for the supplementary catalogue (seelAlexander et al. 2003 



for details). Although this second catalogue is likely to include 
many spurious X-ray sources, it can be used to robustly identify 
fainter X-ray counterparts associated to known sources (see e.g.. 
Sect. 3.4.2. of |Alexander et al.|2003| l. 

These catalogues were used to identify the X-ray counter- 
parts of the 489 VLA/24 /im sources. The 24 fim positions were 
matched to the X-ray positions using a small search radius of 
1.5"; taking into account the high positional resolution of the 
Chandra data (median positional uncertainties »!0.3") and the 
small pixel size of the MIPS-24 fj.m images, the majority of the 
true counterparts are expected to lie within this radius. In fact, 
calculating the probability P to find a random object within 1.5'' 



from the X-ray positions following the prescription of Downes 
|et al.[ ( T986[ ), we obtained a maximum probability of random as- 
sociation P - 0.03 (considering a space density of 24 fim sources 
ofn = 5.2x lO-^deg-^). 

From the main catalogue we found X-ray counterparts (in the 
0.5-8 keV energy band) for 137 of the 489 VLA/24 jum sources 
(a!28%), with a median positional separation of a!0.2". Amongst 
the matched sources, we found that in none of the cases was there 
more than one counterpart within the search radius, with the clos- 
est neighbours being at separations > 2". From the supplemen- 
tary catalogue we identified a further 22 X-ray counterparts to 



Thanks to the large spectroscopic follow-up observations per- 
formed in the GOODS-N field, x!3000 redshift identifications are 
available for the objects in this field. A compilation of spectro- 
scopic redshifts (Zspec) were obtained from the major publicly 
available spectroscopic redshift surveys of the GOODS-N field 
[Chapin et al. 2009; Chap man et al.[ 



(i.e.. 



Barger et al.][2008 
[Cowie et al.|2004| 



Wirth et al.j2004^ , as well as some un- 



published spectroscopic redshift identifications (courtesy of M 
Dickinson). The optical positions of the sources in these cata- 
logues were matched with the 24 fim positions of our detected 
sources using a search radius of 1.0". Since the errors on the op- 
tical positions are typically very small, we used a smaller search 
radius than that used for the X-ray catalogues (see Sect. |2.4| l. 
With this search radius and considering the sky density of the 
24 fj.m detected sources (Sect. |2.4| l, we estimated the spurious 
detections to be ~3%. 

Spectroscopic redshift measurements were found for 1225 
sources amongst the 24 fim detected sample (~63%), with the 
large majority coming from the Barg er et al.[ (2008) spectro- 
scopic redshift catalogue (1030/1225). Two more redshift iden- 
tifications were obtained from Spitzer IRS MIR spectra ( [Murphy 



et al.|2009 i, yielding a total of 1227 Zspec measurements for our 



sources. 



In order to increase the redshift identification completeness 
of the sample, we also included photometric redshifts from a cat- 
alogue built following the procedure described in[Pannella et al. 
( |2009b| l and |Strazzullo et al.| (2010). The photometric redshifts 
(Zphot) were estimated using a PSF-matched multi-wavelength 
catalogue including 10 photometric optical/near-IR passbands 
(from the U band to 4.5 /vm), through a comparison with a li- 
brary of galaxy SED templates, spanning a wide range of galaxy 
types (from elliptical, to star-forming to QSO-dominated) and 
star formation histories (SFHs). The construction of the multi- 
wavelength catalogue and the photometric redshift estimates will 
be described in detail in a paper by Pannella et al. (in prep). 
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Fig. 1: Galaxy stellar masses in units of M© for the entire 24 
fim detected sample with redshift identification (spectroscopic 
or photometric). The stellar masses for the sources detected in 
the VLA radio band are shown as grey histogram and the X-ray 
detected sources are shown as shaded black histogram. 



The photometric redshift catalogue includes 1893 Zpimi within the 
GOODS-//er5c/ie/ area considered here. Photometric redshift 
estimates were available for 67 1 of the 24 fim detected sources 
(~35%) without Zspec measurements. To verify the reliability of 
the photometric redshifts we compared the Zphoi with the spec- 
troscopic redshifts from jBarger et al. ( j2008 ); in the photomet- 
ric catalogue by Pannella et al. there are 1030 sources overlap- 
ping with the Barger et al.| ( 2008) 1 spectroscopic sample. The 
relative accuracy of Zphot, defined as the average absolute scatter 
(AAS = m ean[| Az | /(I + z^pec)], where Az = (Zphoi - ZspecY, e.g. 
[Raffe rty et al.l 2011[) is ^5%, w ith *4% of outliers (AAS > 0.2; 
see also Mullaney et al. 2012 1. Two further Zphot were taken 



from Pop e et al.| ( |2006| l, yielding a total of 673 photometric red- 
shift estimates for our sample. This gives a final redshift iden- 
tification completeness of ~98% (1900/1943 sources, includ- 
ing spectroscopic and photometric redshifts) amongst the 24 //m 
detected sources in the GOODS-Herschel field with a redshift 
range z = 0.02 - 6.54. 

For the purposes of our analysis, we want to investigate 
here only the sources with a significant radio detection, in or- 
der to have reliable measurements of the FIR-radio correlation. 
We therefore only included in our sample sources with a red- 
shift identification amongst the VLA/24 fim sample (484/489 
sources). Due to the limitations dictated by our SED fitting tool 
(Sect. |3.1| l, we also imposed a redshift limit of z < 3.0 to our 
sources, yielding a sample of 458 VLA/24 /im sources with spec- 
troscopic or photometric redshifts of z < 3.0. The analysis of 
the whole 24 yum detected sample, including radio undetected 
sources, will be presented in a future paper (Del Moro et al., in 
preparation). 

2.6. Stellar masses 

The multi-wavelength optical/near-IR catalogue used to estimate 
the photometric redshifts (Pannella et al., in prep.; see Sect. |2.5| l 
was also used to calculate the galaxy stellar masses (M»). The 
stellar masses have been derived using the SED fitting code de- 
tailed in [Drory et aL] ( |2004| |2009| l to fit our multi-wavelength 
data. The star formation histories have been parameterised with 



a linear combination of a main SF event, with SFR exponentially 
declining with time as tp{t) cc exp(-t/T) (where the time-scale 
T = 0.1 - 20 Gyr), and a secondary burst. The main component 
has solar metallicity and an age between 0.01 Gyr and the age 
of the Universe at the source redshift, while the secondary burst 
is limited to <10% of the galaxy total stellar mass and is mod- 
elled as a 100 Myr old constant SFR episode with solar metal- 
licity. We adopted a Salpeter (1955) IMF for both components 
and an extinction law ( jCalzetti et al.| 2000), allowing ranges of 
Ay - 0-1.5 mag and Ay - - 2.0 mag to extinguish the 
main component and the burst, respectively (Pannella et al., in 
preparation; see also |Mullaney et al.|2012[ ). The uncertainties on 
the stellar masses are estimated from the dispersion on the mass- 
to-light ratio {MIL) distribution of the entire library of models 
adopted, as well as from the systematic uncertainties (due to the 
adopted models, IMF, SFH, metallicity, etc.; see e.g. Marchesini 
et al. [200^ 1. The uncertainties are typically larger at low stellar 
masses and range from ^Q.A dex at log M» = 9.0 M© to aiO.2 dex 
at log M, = 1 1.0 Mq (Pannella et al., in prep.). We note that also 
AGN emission in the UV/optical band can cause uncertainties on 
the stellar mass estimates. However, the contamination from the 
AGN affects the stellar mass only when the AGN is very lumi- 



nous (Lx > 10"^^ erg s'^; e.g. [Mullaney et al.|2012| Rovilos & 



|Georgantopoulos|2007tpCueetal.|2010| l 

Stellar masses were measured for 1894 of the 1900 24 //m de- 
tected sources with a redshift identification within the GOODS- 
Herschel field (Sect. |23]l, which include 456 of the 458 VLA/24 
yum sources with spectroscopic or photometric redshift z < 3.0 
that constitute our final sample. The galaxy stellar mass values 
obtained range between log = 7.0 - 12.5 M© (log M, = 
8.7 - 11.9 Mq for the VLA/24 jim sample), with the large ma- 
jority of the sources (a!97%) having log M» =9.0-11.5 Mq 
(Fig. [T]i|^ The median stellar mass for the 24 yum sample is log 
M, Si 10.4 Mq, and for the VLA/24 yum sample the median is log 
Mt » 10.8 Mq. We note that the stellar mass distribution of the 
VLA/24 yum sample sources is consistent with that of the X-ray 
detected sources (median stellar mass: log M, » 10.9 Mq; see 
Fig.[T]). 

3. FIR-radio correlation 

3.1. SED fitting approach 

The emission observed in the MIR and FIR bands is produced 
by dust heated by the radiation emitted through star formation 
and/or accretion onto a SMBH. Star-formation, which occurs 
on large scales in galaxies, heats the dust to a wide range of 
temperatures: the hot dust produces emission at near-IR (NIR; 
A X 2.0 - 5.0 yum) and MIR wavelengths and gives rise to the 
characteristic PAH features (e.g 



Smith 



Cha ry & Elbaz|[200T 

et al.|2007 1, while a large amount of colder dust (Tjust ~ 20 - 50 
K; see also Sect. [TJ produces a typical SED that peaks at FIR 
wavelengths (A ~100 jjm). AGN activity yields on average hot- 
ter dust temperatures (Tdu.it ~ 200 - 1000 K) than star formation, 
so that the bulk of the AGN emission is produced in the MIR 
band with a peak at shorter wavelengths (A ~ 20 - 30 //m; e.g. 
Mullaney et al.|201 l[|Netzer et al.|2007p . The lack of colder dust 
(Tdiist < 200 K) causes a fast decline of the SED at wavelengths 
longer than A > 30 yum (e.g. ,Mullaney et al._2011, Netzer et aL] 
2007| l. AGN activity and star formation are often coupled in a 



^ We note that the fraction of sources with very low (log M, < 8.0 Mq) 
or very high (log M, > 12.0 Mq) mass values amongst the 24 fim de- 
tected sample (3%) is consistent with the estimated fraction of photo- 
metric redshift outliers (see Sect. 2.5 i. 
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Fig. 2: Comparison of the five star-forming galaxy templates 
used in our approach (black lines) with the 105 |Chary & El-| 
|baz| f2001 ) galaxy templates (grey lines); all of the templates 
have been normalised to 1 at 30 yum. With our 5 SFG templates 
we cover a broad range of star-forming galaxy properties, even 
broader than those of the Chary & Elbaz (2001 1 templates. 



galaxy and it is not trivial to separate the emission due to these 
two processes. To disentangle the two components we therefore 
performed a detailed analysis of the IR SEDs of the 458 VLA/24 
fj.m detected sources, with spectroscopic or photometric redshifts 
out to z = 3 (Sect. |23|. 

To represent the galaxy emission we used five star-forming 
galaxy (SFG) templates defined by Mullaney et al. (2011 1, cov- 



ering the wavelength range 6-1000 fim. These five templates are 
defined as composites o f a sample of l ocal star-forming galax- 
ies with LiR < 10'^ Lq ( Brandl et al.| 2006) and are designed to 
sample the full range of IRAS colours observed for these galax- 
ies (see [Mullaney et al.|[2011| for details). We have extended 
the five SFG templates to shorter wavelengths (3 /im) using the 
average starburst SED derived by Dale et al.|( |200T] l. To verify 
whether the |Dale et al.| (|2001) template was suitable for extend- 
ing the Mullaney et al. (2011) templates, we obtained publicly 
available NIR and MIR dat a (from the NASA/IP AC Extragalac- 
tic Database, NEE[^ for the [Brandl et al.| ( |2006| l sample of local 
star-forming galaxies and we plotted these data points over our 
extended templates as a check; we thus verified that data and 
models agreed with reasonable scatter. 

We also extended the five SFG templates to the radio band 
with a power-law slope Sy oc v with a = 0.7 (e.g., |Ibar et al. 
|2010|[2009 ; see Fig. [2]). The normalisation of the radio power- 
law component was fixed according to the typical FIR-radio rela- 
tion found for local star-forming galaxies ( |Condon|1992[ |Helou| 
|et al.|[T985] l. We stress that even though we limited our SFG 
template library to only 5 templates, they represent a wide range 
of IR color-color properties of star-forming galaxies ( |Mullaney| 
et al.||2011,), ev en broader than those reproduced by, e.g., the 
Chary & Elbaz (2001 1 galaxy template library (105 templates; 



To reproduce the emission from the AGN we used the em- 

com- 



pirically defined AGN template by Mullaney et al. (2011 



posed by a broken power-law, declining at wavelengths longward 
of /t > 30 yum as a modified black-body. We note that this tem- 
plate is in agreement with the typical SEDs produced by clumpy 

^ http://ned.ipac.caltech.edu/ 



torus models (e.g. Nenkova et al.|2008a|b I. We fixed the power- 
law indices at the average values Fi = 1.7 and r2 = 0.7, with 
a break point at Agrt - 19 fJ.m (see [Mullaney et ar]|2011| l. We 
have allowed the AGN component to be modified due to dust ex- 
tinction, using the extinction law of Draine ( 2003) , which mainly 
affects the template at /I < 30 fim and also produces the typical 
silicate absorption feature at 9.7 jum, often observed in AGN (e.g. 
Martm ez-Sansigre et al. 2008 ; Ro che et al.|19 91, 2007; ShfetaL 
I2OO6I1. 



The flux densities at 8 fim, 16 fim, 24 //m from Spitzer and 
100 jiva, 160 fim. and 250 //m from Herschel have been used in 
the SED fitting process to constrain the SEDs of our sources (Fig. 
[3]). In the case of non-detections (S/N<3), the measured fluxes 
at each source positions, with the large associated uncertainties, 
were used in the SED fits. The fluxes in the shorter Spitzer- 
IRAC bands (3.6, 4.5, 5.8 /im) were not included in the SED 
fitting process as they fall out of the wavelength range covered 
by our templates at relatively low redshifts z > 0.2. Moreover, 
at these wavelengths the observed SED is often dominated by 
the old stellar population emission, which is not accounted for in 
our SED templates. Data at longer wavelengths from Herschel- 
SPIRE 350 ^im and 500 yum (Elbaz et al. 201 1 ) were also not in- 
cluded in the fits, because of the larger uncertainties on the mea- 
sured fluxes due to strong blending (see Sect. |2.2| l and the low 
number of significantly detected sources. However, these fluxes 
(or upper limits) have been included when plotting the SEDs, as 
a visual check of the best-fitting SED solutions (see Fig. [3] and 



Fig. B.l 



We note that only using photometry at A < 250 pm (SPIRE- 
250 yum band) means that at high redshifts iz > 1.3) we are not 
able to fully constrain the FIR SED peak at wavelengths longer 
than /I X! 100 pm (rest-frame) and therefore we cannot exclude 
the presence of a colder dust component, which would yield 
higher FIR fluxes than that predicted from our model. However 
we can anticipate (see Sect. [3.2[ ) that we find very good agree- 
ment between our average FIR-radio relation and that found in 
several previous works. 

The SED fitting process was developed as follows: 

1. Initially, only the star-forming galaxy templates were used in 
the fit; we fitted the data of each source with each of the five 
SFG templates (Fig. |2j using minimisation to evaluate the 
goodness of the fit; 

2. As a second step, we performed new fits, again using 
minimisation, by adding the AGN component, also includ- 
ing extinction of the AGN component as a free parame- 
ter (varying between Ay ~ 0-30 mag, corresponding to 
N}\ a: - 5 X 10^~ cm"^ assuming the av erage galactic dust- 



to-gas ratio Ay = A^h/(1-8 x 10^'); e.g. [Predehl & Schmitt 
T9M1), to each SFG template (SFG + AGN); 



3. Finally, an /-test was performed using the;i'^ values and the 
degrees of freedom (d.o.f.) for all of the five pairs of solu- 
tions to evaluate the improvement of the fit due to the addi- 
tion of the AGN component. We accepted the SFG + AGN 
model as the best-fit if the AGN component significantly (i.e. 
> 90% confidence level, according to the /-test probability) 
improved the resulting for the majority of the solutions, 
i.e. in at least three out of five fitting solutions. 

The criteria we adopted to define the best-fitting model were es- 
tablished after performing several tests on our SED fitting ap- 
proach; in particular, we tested these criteria on a sub-sample 
of our sources for which detailed Spitzer-IRS MIR spectra are 
available (see Appendix A, for details). 
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(Sect. 



Examples of spectral energy distributions (SEDs) to demonstrate the variety of SEDs found for the radio-excess sources 
The SEDs on the left are consistent with a simple galaxy template: SFG template ("IR SFG" 
bottom; 



3.2i 



top) or elliptical galaxy 
We note that the elliptical 



template ("passive", bottom; see Sect. |4.3| i; the SFG template upper limit (grey line) is also shown 
template (long dashed line) is not fitted to the data, but it is only shown to demonstrate that it can well represent the data. The SEDs 
on the right are fitted with a star-forming galaxy (dashed line) -i- AGN (dotted line) model ("IR AGN"): on the top there is an X-ray 
detected AGN, on the bottom an X-ray undetected AGN. The total SEDs are shown as black solid lines. The filled circles represent 
the Spitzer 8, 16, 24 //m and the Herschel 100, 160, 250 yum flux densities (in mJy), used to constrain the SEDs. The open symbols 
indicates the data that were not included in the SED fitting process: red triangles are Spitzer-l^AC 3.6, 4.5, 5.8 /im flux densities, 
black circles are SPIRE 350 and 500 /im and black squares are VLA 1.4 GHz flux densities; the radio data do not match the SEDs 
in these cases, because the sources have excess radio emission compared to that expected from pure star formation (Sect. |3.2| l. The 
blue star represents the 6 luminosity of the AGN predicted from the X-ray luminosity in the rest-frame 2-10 keV, using the |Lutz 
|et al.| ( [2004| relation for local unobscured AGN; we note that this point does not match the IR AGN component because the X-ray 
luminosity tends to underestimate the intrinsic AGN power if the AGN emission is heavily absorbed (see Sect. |5.1| l. On the top left 
corner of each plot the source redshifts are indicated as well as the corresponding source number in Table 1 (col.l). 



Once the best-fit model was defined, the final measurements 
of the properties of the sources (i.e., FIR flux, relative AGN/SFG 
contribution, etc.) and their errors were derived as weighted aver- 
ages of the values obtained from the five best-fit model solutions 
(see Appendix A for details). For the AGN -i- SFG models we in- 
cluded in the average only the solutions where the AGN compo- 
nent was significant (> 90% confidence level). This is because, 
due to the sparse data used to constrain the SEDs, in some cases 
the five solutions obtained from the different SFG templates were 
similar (small difference in;^'- values) and did not allow us to un- 
ambiguously determine a unique solution that best characterised 
the data (see Appendix A, for details). 



3.2. Radio-excess sample selection 

The FIR-radio correlation is typically defined as the ratio be- 
tween the flux in the FIR band (~ 40 - 120 |/m; rest-frame) and 



the flux density in the rest-frame radio band ( 1 .4 GHz; e.g. Helou 
|et al.|1985||Sargent et al.|2010 ai Many studies performed so far 
at high redshift relied on either MIR fluxes as a "proxy" for the 
FIR (or bolometric IR) flux (e.g. from the S24 um', Appleton et al. 
IDonley et al.||2005j, 70 fj.m flux density (e.g. ,Appleton 



2004 



et al. 



et al. 



2004; B ourne et al.|201 1| [Sargent et al.|[20TDal |Seymour| 
2009 j l, or on SED fitting spanning only the rest-frame MIR 



band (e.g. [Padovani et al. 2011 Sargent et al.||2010a| . These 
methods are often inaccurate because they require assumptions 
about the source SED over the whole IR band and/or on the bolo- 
metric corrections. Through our detailed SED analysis of the 
458 VLA/24 fim sources with z < 3.0, using the approach de- 
scribed in the previous section (Sect. |3.1| l and the Herschel data 
to constrain the FIR SED peak, we can overcome these issues by 
directly measuring the galaxy emission over the whole IR band. 

From the best-fit models, we calculated the FIR flux (/fir) 
by directly integrating the total SEDs over the rest-frame wave- 
length range A = 42.5 - 122.5 jum (e.g. [Helou et al.7l985 ). We 
used the total SED, which in many cases includes contributions 
from both SFG and AGN components, to calculate the FIR flux 
because we aim at a conservative selection of radio-excess AGN, 
since radio quiet AGN often follow the typical FIR-radio relation 
of star-forming galaxies. 

In only ^3% of the 458 analysed sources (15 sources) was 
the SED fitting analysis unable to provide a good representation 
of the mid- and far-IR dataj^In most of the cases (6/15), this was 
due to large uncertainties on the photometric redshift, or pos- 
sibly to spurious counterpart associations between the different 



The poor fits were flagged by visual inspection of the resulting SEDs; 
these are the cases where even the best-fitting SEDs deviate significantly 
from the data. 
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Fig. 4: Left: Distribution of FIR-radio flux ratio {q) for afl of the sources in the VLA/24 jum sample (458 sources, excluding the 
upper limits; see Sect. |3.2| i; the blue Gaussian profile represents the best fit to the peak {qmean - 2.21 + 0.18), while the red dashed 
Gaussian is the q distribution found by |Helou et al.| ( [T985| l for local starburst galaxies. The vertical fine indicates the separation 
we assumed to select our radio-excess sources, q - 1.68, corresponding to 3cr from the peak (2cr when include radio upper limits). 
Right: FIR-radio flux ratio (q) as a function of redshift for all of the VLA/24 yum sources; the radio-excess sources have q < 1.68 
(dot-dashed line); the solid horizontal line represents the average q value for the "radio-normal" sources {q - 2.21) and the shaded 
region indicates the +lcr error from this average. The radio-excess sources identified by Donley et al. ( 2005| l are plotted as crosses 
(see Sect. |3.3| l. The binned q average for the "radio-normal" sources, in three redshift bins, is also shown as black circles; we found 
no significant evidence of evolution of the FIR-radio correlation out to z = 3. 



catalogues (see Sect. 2.5 i. In a smaller number of cases (4/15), 
the poor fitting results seemed to be due to large uncertainties on 
the flux density measurements, especially at wavelengths where 
the sources are not significantly detected (S/N < 3). We flagged 
these "problematic" cases in Table 1 (column "Fit"), being aware 
that the measurements obtained from their SEDs are not fully re- 
liable. The remainder of these sources (5/15) have strong emis- 
sion in the Spitzer IRAC bands, even stronger than the flux de- 
tected at 24 jum (Sii/S 24 > 1), and are not detected at longer 
wavelengths (FIR) by Herschel. This suggests that the emission 
from star formation (or AGN) in these sources is weak, while 
the emission from old stellar population, usually dominating the 
rest-frame NIR band (/I < 5 //m, rest-frame), is strong; their in- 
frared SEDs are therefore more consistent with that of galaxies 
dominated by passive stellar populations rather than by active 
star formation. For these sources FIR flux upper limits were esti- 
mated by normalising the SFG templates to the 24 jim datapoint 
(Fig. [3]and Fig. |B.1[ ), which is thus used as a proxy of star forma- 
tion, and integrating the SED between 42.5 - 122.5 ;um. Of the 
five measurements obtained from the different SFG templates, 
the maximum has been taken as the /fir upper hmit. 

Using the definition given by |Helou et al.| ( [T985| l, we calcu- 
lated the ratio between the far- infrared and radio emission (q) as: 



q = log [/fir/(3.75 x lO'^ Hz)] - log [5,(1.4 GHz)] 



(1) 



where /fir is in units of W m"^, 3.75 x 10'^ Hz is the frequency 
at the centre of the FIR band {A - 80 jum) and 5,(1.4 GHz) 
is the radio flux density (in units of W mT^ Hz"') at rest-frame 
1.4 GHz, extrapolated from the VLA data using the power-law 
slope Sy oc v~", with a = 0.7, typical for star-forming galaxies 
(e.g., |Ibaret al.|2010| [20091 ) . 



The q distribution obtained for the entire sample is shown in 
Fig. [4] (left); the sources where /fir (and therefore also q) is an 
upper limit are not included in the histogram. The peak of the 
distribution is at ^ a; 2.2, in excellent agreement with the values 
typically obtained for star-forming/starburst galaxies (e.g. Coi>] 
|don||1992t [Helou et aL]|1985| l. However, the distiibution is not 
symmetrical around the peak and shows a broad tail at low q val- 
ues, indicating a relatively large number of sources with excess 
radio emission over that expected from star formation processes. 

In order to identify the radio-excess sources, we fitted the 
peak of the q distribution, considering only sources with q > 2.0, 
with a Gaussian profile and estimated the spread of the FIR-radio 
correlation expected for star-forming galaxies: from the best-fit 
we obtained q - 2.21 ±0.18 (see Fig. [4]i. If we also include the 
radio undetected sources (S/N < 3) in the q distribution, we ob- 
tain a Gaussian profile with a very similar peak, but larger scatter 
(q - 2.24 + 0.29), due to the larger uncertainties on the q values. 
The q values for the radio undetected sources were calculated us- 
ing the radio 1 .4 GHz flux measurements at the 24 yum source po- 
sitions (Sect. [2.3| l. We set the separation between "radio-normal" 
and "radio-excess" sources at ^ = 1.68, corresponding to a 3cr 
deviation from the peak of the distribution for the VLA/24 //m 
sample (~2cr from the peak for the whole 24 fim sample). We 
defined "radio-normal" as the sources with q > 1.68, that follow 
the typical FIR-radio relation {q x 2.2); we note that this popu- 
lation includes star-forming galaxies, but also typical radio-quiet 
AGN ( e.g. [Moric etal.|2010[ [Padovani et al.|2011[ [Sargent et aL 



[2010a[ ), while the large majority of the radio-excess sources are 
most likely to host AGN activity. 

With our selection we obtained a sample of 5 1 radio-excess 
AGN candidates, ~ 1 1 % of the whole VLA/24 yum detected sam- 
ple at z < 3.0. In Fig. [4[ (right) we show the FIR-radio ratio 



Article number, page 8 of 



28 



A. Del Moro et al.: Radio-excess sources in GOODS-Herschel. 



2 - 







-1 - 



-2 

3 



° 1 



-1 



% ® 



® 



® 



® 



® 



I I ' I ' I 




• Radio detected 

® Our radio-excess (q<1 .68) 



0.0 0.5 



1.0 



1.5 2.0 
Redshift 



2.5 3.0 



Fig. 5: Infrared-to-radio flux ratios, using the observed 24 jim 
iq24', top) and 100 /im flux densities (^loo; bottom); the radio- 
excess sources in our sample < 1 .68) are highlighted by open 
circles. The shaded regions in the two plots represent the ^24 
and ^100 ratios predicted for our range of star-forming galaxies 
as a function of redshift. The horizontal hne in the top plot is 
the q24 - selection used by| Donley et al.| ( [2005| l to define their 
radio-excess sample. In the bottom panel we identified ^100 = 1.5 
as an easy selection of radio-excess sources (dashed line), with 
reasonable completeness (=s60-80%) at least out to z ~ 2. 



(q) for the entire sample as a function of redshift. The horizontal 
dashed line indicates the q value at the peak of the distribution 
(q = 2.21 ± 0.18) with Icr uncertainty (shaded region), while the 
dot-dashed line at ^ = 1 .68 represents our threshold for radio- 
excess sources (3cr deviation from the peak). In the plot we also 
marked 12 of the IDonley et al. (2005 1 radio-excess sources that 
are found within the GOODS^Herschel (North) field and overlap 
with our sample|^(see Sect. 3.3 1. In Fig. |4](right) we also show 
that the average q for all of the radio-normal sources, calculated 
in three diff'erent redshift bins (zi - 0.0 - 1.0, Z2 = 1-0 - 2.0, 
Z3 - 2.0 - 3.0, black circles), remains fairly constant, within the 
errors, over the whole redshift range, and therefore, the apparent 
decrease of the q values at high redshift is not significant (e.g. 



' We note that for 6 sources the redshifts listed in Donley et al. 1 2005 1 
are in disagreement with ours (see Table [TJ: in five cases we have new 
spectroscopic redshifts where there were only photometric redshift esti- 
mates (or no estimates at all) in |Donley et al.| ( (2005| l; in the remaining 
case, we have a photometric redshift where there was no redshift mea- 
surement for this source in [Donley et al.]j2005} . 



Bourne et al.||20lTHEibaz et al.|2002[ |Mao et al.|[20n) [Sargent 



et al.|2010bl l. 



3.3. IR-radio relation: 

q24 and q\QQ 

Several previous works that have investigated the FIR-radio cor- 
relation and radio-excess sources have used different methods to 
define the FIR-radio flux ratio. In particular, in many studies the 
FIR flux has been replaced by the monochromatic flux density at 
24 yL(m (e. g. [Appleton et al |2004[ [Donley et al.|2005t[Ibar et a 



[2008; Sarg ent et al. [[2010a | l , 60 urn (e.g. JVlahakis et al.[[2007 



or 70 yum (Mao et al. 201 1 ; S argent et al.|201()a Seymour et a 
[2009 ), as a proxy for the FIR emission of the galaxy. It is there- 
fore interesting to see how these definitions of the FIR-radio ra- 
tios, calculated from monochromatic flux ratios, compare to q es- 
timated by us (Eq. 1) across the full FIR waveband. In particular 
we performed a direct comparison of our sample selection with 
that used by Donley et al.|(200 5), where the radio-excess sources 
were selected using qiA < 0, with ^24 = log (^24 1.4 ghz), 
calculated using the observed 24 /vm and 1.4 GHz flux densities. 
As shown in Fig. |4] (right), the majority of the [Donley et al 



(2005) sources (11/12 that overlap with our sample) are also 
radio-excess sources based on our definition. However g'24 does 
not provide a complete selection, since the majority of our radio- 
excess sources (a!78%, 40 objects) are missing from the [Donley] 
[et al.| (2005) sample. To verify that this is not only due to the 



improved sensitivity of our radio data (Sect. 2.3 1, we compared 
the number of radio-excess sources we would identify with our q 
definition and selection (^ < 1 .68) using the radio catalogue by 
'Richards (2000), the same as that used by [Donley et al. (2005^: 
we obtained 24 radio-excess sources in the GOODS-Herschel 
field, twice the number of those identified by |Donley et al.[ ( |2005| l 
using ^24 within the same sky area (12 sources). Moreover, 
in one case the q24 selection by [Donley et al. ( 2005 [ ) disagrees 
with our FIR-radio ratio, underestimating the FIR power of the 
sources (Fig. |4] right). This effect is shown more clearly in Fig. 
|5] where we calculated the observed g'24 and qioo for all of the 
VLA/24 yum detected source^ The q24 and ^100 ratios predicted 
for our range of star-forming galaxies are shown as shaded re- 
gions in Figure|5] where our radio-excess sources (q < 1.68) are 
marked with open circles; the ^24 = line adopted by [Donley] 
et al. (2005 1 to select radio-excess sources is also shown (Fig. ]5J 
top panel). Although the ^24 < criterion would select some 
of our radio-excess sources, more than 60% of them would be 
missed, testifying the large incompleteness of this method when 
compared to the FIR-radio relation q (Eq. 1). Moreover, using 
q24 < to define the radio-excess sample would include some 
"radio-normal" galaxies (according to our definition; Sect 3.2 1, 
introducing some contamination to the sample (~10%). 

On the other hand, the q'loo ratio (Fig. ]5] bottom) seems to 
be in better agreement with our q selection than the ^24 ratio, as 
few of our radio-excess sources fall into the shaded area. This is 
expected since the flux observed at 100 fim is closer to the FIR 
SED peak than that observed at 24 fim which is also likely to be 
affected by the PAH emission features at z > 1 (as illustrated by 
the shaded region in Fig. ]5] top) and/or by AGN contribution. 
By defining a separation line below the shaded area to identify 



* qioo = log (S 100 ^m/S,A GHz), where S m and S 1.4 ghz are observed 
monochromatic fluxes. Since we do not have photometric measurements 
at 60 fim or 70 fim to directly compare our results with those obtained 
in previous works (e.g. [Mao et al.|201 l[[Sargent et al.|2010a[[Seymour[ 
[et al.|2009[ [Vlahakis et al.|2007| )' we used the observed 100 fim flux 
density, which is the closest available data point to allow any comparison 
with qeo or qjo. 
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Fig. 6: Top panel: Fraction of radio-excess sources as a func- 
tion of the rest-frame FIR luminosity (42.5-122.5 //m, in units 
of L0); the fraction decreases with increasing Lfir, as expected 
since sources with high Lfir are more likely to be luminous star- 
forming galaxies. Bottom panel: Rest-frame FIR vs. radio 1.4 
GHz luminosity (rest-frame); the radio-excess sources are plot- 
ted as red circles and some local well known Compton-thick 
AGNs (open squares) are shown for comparison (data taken from 
literature). The black dotted line correspond to (^rl) = -0.38 



(average for a sample of RL AGN from Evans et al.( 2Q05); the 
black solid line correspond to ^ = 2.21 (average for "radio- 
normal" sources) and the red dashed line correspond to <7 = 1.68, 
our selection limit for radio-excess sources. 



radio-excess sources (e.g. g'loo = 1-5; Fig. |5]l we would recover 
the majority of our radio-excess sources (~60%; ~80% includ- 
ing the upper limits that lie right above the separation line), a 
much larger fraction than those identified using qiA- We stress, 
however, that the ^100 (observed) selection is still not complete 
and much less accurate than our selection method using the rest- 
frame FIR flux to calculate q. 



4. Results 

4.1. Sample properties 

We established in Section |3] that the direct measurement of the 
FIR flux from the best-fit SED provides a reliable and unam- 
biguous way to identify sources that have radio emission in ex- 
cess to that expected purely from star formation. As described in 
the previous sections, from our analysis we identified 51 radio- 
excess sources, which are most likely hosting an AGN. 

In Fig. |6] we show the rest-frame FIR and the radio 1 .4 GHz 
luminosities (Lfir - 4 ghz; rest-frame) for all of the "radio- 
normal" (black crosses) and the radio-excess sources (red cir- 
cles) in the GOODS-//er5c/ie/ field. The diagonal lines repre- 
sent the mean q value obtained for the "radio-normal" sources 
(q = 2.21) and the separation between "radio-normal" and radio- 



excess sources (q - 1.68). For comparison we also indicate 
in the plot the average q obtained for a sample of low-redshift 
(z < 0.15) radio-loud AGN (RL; <^rl) = -0.38) from [Evans 
|et al.| f2005 ). For these sources the FIR flux was calculated us 



ing the photometry at 60 /urn and 100 //m, following Helou et al. 
( |1985| l. So me well studied local C ompton-thick AGN taken from 



literatur^( |Della Ceca et al. 2008 1 are also plotted, as a guide; al- 
though the majority of these sources have FIR-radio ratios con- 
sistent with that of star-forming galaxies and most radio-quiet 
AGN (q K 2.2), three out of the eight AGN are radio-excess 
sources (i.e., Mkn 3, NGC 7674 and NGC 1068), suggesting that 
some Compton-thick AGN might also be included in our radio- 
excess sample. 

It is important to note that, although some of our radio- 
excess sources have radio luminosities typical of radio-loud 
AGN {Lrad > 10^"* W Hz"'), the majority lies in a region in be- 
tween those occupied by RL AGN and "radio-normal" sources 
and therefore can also be referred to as "radio-intermediate" 
sources (e.g. [Drake et al.|2003| l. We therefore stress that the def- 
inition of "radio-excess" does not necessarily mean "radio-loud" 
(see Sect.[r|l. 

In the top panel of Fig. [6j the fraction of radio-excess sources 
as a function of the rest-frame FIR luminosity is shown (with 
Icr uncertainties; [Gehrels|[1986| l. The sources in our sample 
span a wide range of FIR luminosities, Lfir ~ 10^ - 10'^ L© 
(Table [TJ, typical of normal star-forming galaxies, luminous IR 
galaxies (LIRG; Lfir ~ 10" Lq) and ultra-luminous IR galax- 
ies (ULIRGs; Lfir ~ lO'^ Lq). However, the fraction of radio- 
excess sources decreases with increasing FIR luminosity, and 
only one source out of 51 reaches the high luminosities typi- 
cal of ULIRGs. This is not unexpected, since at the high FIR 
luminosities of ULIRGs, the radio core must be very bright 
(Liad > 2 X 10^'* W Hz '; RL AGN regime) to be identified as 
radio-excess above the very strong star-formation activity. 

4.2. X-ray and radio emission from the AGN 

Since the excess radio emission is attributed to the presence of an 
AGN, it is interesting to investigate the X-ray properties of these 
sources. Amongst our radio-excess AGN sample, only ~53% 
(27/51) are detected in the X-ray band (Table 1), hence about 
half of the sample is X-ray undetected. This suggests that the 
radio-excess selection is a powerful method to identify AGN that 
would otherwise be missed by even the deepest X-ray AGN sur- 
veys, such as the Chandra 2 Ms X-ray survey ( [Alexander et aL] 
[2UU3] l. 

In Figure [7j we show the fraction of radio-excess sources 
as a function of X-ray luminosity. The sources have been di- 
vided into five X-ray luminosity bins from L2_io kev - lO''" erg 
s"' to L2-iokeV - 10"*^ erg s"'. We found that the fraction of 
radio-excess AGN increases from ^7% at low X-ray luminosi- 
ties (L2_io keV < 10"*^ erg s"') to ; 
bin(L2-iokeV = 10^4 -1045 ergs- 
excess AGN at high luminosities suggests that the brightest AGN 
also produce the most powerful radio emission (e.g. Brin kmann[ 
et al.[[200"0 La Franca et al.[[2010j ), presumably the majority of 



k60% in the highest luminosity 
-'). This large increase of radio- 



it coming from a radio core. We note that the fact that more lu- 
minous AGN are also more radio bright does not mean that high 
luminosity AGN are typically radio loud. On the contrary, the 
radio-loudness, measured as Rx - log (vLi.4 ghz/L2~io kev), de- 
creases with increasing X-ray luminosity (e.g. [La Franca et al.[ 
2010). 



' The FIR and radio data for these sources are taken from the 
NASA/IPAC Extragalactic Database (NED). 
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Fig. 7: Fraction of radio-excess sources in different X-ray lumi- 
nosity bins; the fraction of sources with excess radio emission 
increases with X-ray luminosity, from ^1% at Lx - 10""' - IC*' 
erg s"' to -60% at Lx = 10^4 io45 g^g s" '. The eiTor bai's 
correspond to Icr uncertainties ( Gehrels|1986[ l. The dashed his- 
togram represents the radio-excess fraction including the X-ray 
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Fig. 8: Radio core flux density at 1.4 GHz measured from the 
VLBI data (Chi et al. 2012, submitted) vs. the excess radio 
emission (i.e., S'i^ghz, tot - 1.4 ghz, sfg) for our radio-excess 
sources. Our radio fluxes agree with the VLBI measurements 
within a factor x 2; our fluxes are typically higher because they 
might include some extended radio emission from outflows, or 
"aborted jets" (e.g.|Ghisemm et al.|20"04|. 



To investigate in detail the radio properties of our radio- 
excess sources we compared the excess radio emission (i.e. the 
excess above the radio emission expected from star formation), 
with the AGN radio core emission detected in deep VLBI 1.4 
GHz obsei-vations of the Hubble Deep Field North (HDF-N) and 
the Hubble Flanking Fi elds (HFF) The details of the VLBI ob- 
servations are given in |Chi et~aL ( 2009) and Chi et al. (2012, 
submitted). Briefly, the sensitivity of the VLBI data varies sig- 
nificantly from the centre of the field (rm.s. noise level of 7.3 
fiiy/heam within 2' from the phase center) to the outer parts 
of the field (2' - 8'; rm.s. noise level of 14—37 yuJy/beam), 
which means that only sources with a radio core brighter than 
S 1.4 GHz ^100 yuJy are likely to be detected. Amongst our radio- 
excess sample, only 13 sources are bright enough (on the basis 
of the VLA 1.4 GHz flux density; Table [1]) to be detected in the 
VLBI images and for eight of them a compact radio core was in- 
deed detected by VLBI, positively confirming that these sources 
host an AGN and that the excess radio emission is in fact due to a 
radio core. The remaining five sources have on average lower ra- 
dio VLA fluxes compared to the sources detected by VLBI (the 
only exception is #47, in Table [1] which is an extended radio- 
loud, wide-angle-tailed source) and lie in the outer regions of the 
VLBI field; the non-detection therefore might be due to sensitiv- 
ity issues, or, as in the case of source #47, to the radio flux being 
dominated by extended emission, more than a compact core. The 
rest of our radio-excess sources lie outside the VLBI region (13 
sources) or are typically too faint (S 1.4 ghz < 100 fiJy) to be sig- 
nificantly detected above the VLBI sensitivity limit (25 sources). 



In Figure [8] we show the comparison between the AGN ra- 
dio core emission measured from the VLBI data and the excess 
radio emission obtained for our sources. The excess radio emis- 
sion was estimated by subtracting the 1.4 GHz flux predicted for 
star formation assuming the average value q - 2.21 (Sect. 3.2 1, 
from the total radio flux density (from the VLA data at 1.4 GHz; 
Sect |2.3| l. The obtained radio excess flux densities agree with the 
VLBI AGN core emission within a factor of x2. This is a further 



validation that our SED analysis is reliable and that with it we 
are able to predict the AGN radio-core luminosities with reason- 
able accuracy. We note that our estimates tend to slightly over 
predict the AGN radio-core flux densities because they might in- 
clude contributions from extended radio emission, such as com- 
pact lobes, or "aborted jets" (e.g. GhiselUni et al.|2(j04 l, which 
are invisible for the high resolution VLBI. 



4.3. SEDs of the radio-excess AGN 

From the SED fitting analysis we found that the infrared SEDs 
of the radio-excess sources are rather varied. A significant 
AGN component, which typically dominates the emission below 
~40 i-im (MIR), has been detected in ~45% of the radio-excess 

this fraction is 



sources (23/51; hereafter, IR AGN; see Fig. B.l 



much higher than that found for the "radio-normal" population 
(~8%; Del Moro et al., in preparation), in line with the fact that 
radio-excess sources are AGNs, as opposed to "radio-normal" 
sources, which are a mix of radio-quiet AGN and star-forming 
galaxies (Sect. |3.2[ ). The average FIR-radio ratio for the IR AGN 
amongst the radio-excess sources is (^)ir agn - 1.02 + 0.11, 
significantly lower than our threshold adopted to separate radio- 
excess from radio-normal sources (q < 1.68). The remainder of 
the sample, however, do not require a significant AGN compo- 
nent in the SED fitting: ai45% (23 out of 51 sources; Table |2| 
have SEDs consistent with those of star-forming galaxies (here- 
after, IR SFGs; see Table [T] and Fig. |B.l| i and the remaining 
~10% (5 out of 51) have IR SEDs not compatible with either star 
formation or AGN activity (see Sect. |3.2| l. For the latter group 
of sources we verified that the IR data are better represented by 
a normal elliptical galaxy template (generated with the GRASIL 
code; |Silva et al. |1998|l, suggesting that these sources might be 
"passive" (not star-forming) galaxies (see Fig. B.l 1. The mean q 
value for the IR SFGs is <^)ir sfg = 1.33 + 0.09, larger than the 
average obtained for the IR AGN, but still well below the radio- 
excess selection limit {q < 1.68). For the "passive" sources, the 
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Table 2: Summary of the IR and X-ray properties of the radio-excess sources. 



Class 


Total 


% 


X-ray detected 


% 


X-ray undetected 


% 


{q) 


(SFR) (Mo yr-') 


Radio Excess 


51 


100% 


27 


53% 


24 


47% 


1.10±0.08 


58.8+13.8 


IR AGN 


23 


45% 


18 


78% 


5 


22% 


1.02+0.11 


45.5+14.4 


IR SFG 


23 


45% 


7 


30% 


16 


70% 


1.33+0.09 


69.2+21.7 


Passive 


5 


10% 


2 


40% 


3 


60% 







FIR-radio ratio calculated from the SEDs are all upper limits and 
therefore we cannot provide a reliable average q value. 

In the far-infrared band, above A ^ AQ /vm, the emission 
is dominated by star formation peaking around 100 /im (rest- 
frame) for the majority of the radio-excess sources (a;86%; 44 
sources); this fraction also includes the majority of the sources 
that are dominated by AGN emission at shorter wavelength (i.e., 
IR AGN). However, we found a number of sources (^14%; 7 
sources) where the AGN emission outshines the star formation 
over the whole IR band, with the SFG contributing less than 50% 
to the total FIR luminosity; these sources are amongst the bright- 
est AGN in our radio-excess sample and the majority of them 
(6/7) are also detected in X-rays. Constraining the actual contri- 
bution from star formation (if there is any) to the total SED in 
these objects would require deep sub-mm observations, in order 
to sample the wavelengths beyond the rest-frame 100 jum. 

Using the broad-band IR luminosities (/I = 8 - 1000 /im, rest- 
frame) of the star formation component (i.e. removing the AGN 
contamination) derived from the best-fit SEDs in the same way 
as for the FIR luminosities (Sect. |3.2[ ), we estimated the star for- 
mation rates (SFRs) for our radio-excess sources. To convert the 
IR luminosity (Lir) into SFR we used the relation from Kenni- 
|cutt| ( |1998[ ), assuming a Salpeter (1955) IMF. The SFRs obtained 
for our sample span a wide range of values, from SFR ^ 0.8 
Mq yr"', indicating very low star formation, to SFR » 350 M© 
yr~', typical of starbursting systems. The average star forma- 
tion rate of our radio-excess sources is (SFR) = 58.8 + 13.8, 
with the IR SFG having slightly higher values than the IR AGN 
«SFR)iR SFG = 69.2 + 21.7 and <SFR)ir agn = 45.5 + 14.4; Ta- 
ble [2]). We note that the SFR could not be calculated for four 
of the IR AGN where the AGN component significantly dom- 
inates the total IR SED (>90% contribution; i.e. #8, #13, #20 
and #21, see Fig. B.l i and therefore the star formation compo- 
nent was not constrained from the SED fitting. Moreover, the 
"passive" sources have not been included in the SFR average for 
the radio-excess, since the SFRs estimated for these sources are 
upper limits. 

We stress that in our calculation of the SFRs we removed any 
contribution from the AGN to the IR luminosity, thanks to our 
SED decomposition. It is important to note that if we did not 
acknowledge the fact that AGN can contribute significantly also 
in the FIR band, and therefore assume that the total IR luminosity 
is only due to star formation (see e.g., Elbaz etal.|2 010; Mullaney 
let al.|2012[ [Santini et al.|2012[ |Shao et al.|20T0l i; we would over 
predict the SFR up to a factor of ~2 when the AGN emission 
contributes <50% to the total SED, and up to a factor of ~4, 
when the AGN dominates (>50% contribution) the SED at FIR 
wavelengths. 

5. Discussion 

We have shown in the previous sections that through detailed 
SED analysis in the IR band using deep Spitzer and Herschel data 
we were able to identify a relatively large fraction of radio-excess 



AGN from an initial radio-IR selected sample in the GOODS- 
North field. Although the excess radio emission is attributed 
to an AGN, we found that the multi-wavelength properties of 
these sources are rather heterogeneous; from the results obtained 
for the radio-excess source SEDs (Sect. |4.3[ ), we can distin- 
guish three main types of objects amongst our sample: i) infrared 
AGN, ii) infrared star-forming galaxies (IR SFGs) and iii) "pas- 
sive" systems. We aim here to explore in detail the variety of 
sources that compose the radio-excess population and whether 
the radio-excess AGN are different from the general X-ray se- 
lected AGN population. 

5.1. Infrared AGN 

About 45% of our radio-excess sample is made of sources that 
we can define as "typical AGN": these sources show a clear AGN 
component at IR wavelengths, which contributes more than 50% 
to the total emission in the MIR band, and are generally also 
detected in the X-rays (a!78%; Table[2|i: they are typically X-ray 
bright AGN, with luminosity above Lx > lO"*"* erg s"'. 

In Fig. [9] (top panel) we show the fraction of IR detected 
AGN at diff'erent X-ray luminosity bins. Although the number 
of sources in each bin is relatively small, it is evident that the 
fraction of AGN detected at IR wavelengths decreases with de- 
creasing luminosity of the sources. While at Lx > 10'*"' erg s"' 
we are able to identify almost all of the X-ray detected AGN 
(>90%) through the IR SED analysis, at lower X-ray luminosi- 
ties (Lx = 10"*^ - lO'*^ erg s"') this fraction drops to a;50%. This 
means that when the AGN is intrinsically bright its emission is 
strong enough to outshine that of the host galaxy and it can be 
clearly detected at IR wavelengths (at least in the MIR band); 
for the brightest sources the AGN dominates the whole IR band, 
even at far-infrared wavelengths (e.g. #8 and #13, see Fig 
Sect.[43|i. 

From the best-fit SEDs we can calculate the MIR luminos- 
ity al A - 6 fj.m for the AGN (extracted from the detected AGN 
component only), which is often used as a diagnostic for the in- 
trinsic pow er of the AGN (e.g. [Alexander et al.||200 8', 'Gandhjj 
|et al.|2009[|Georgantopoulos et ail.|201 l[|Lutz et al.|20 04), since 
the IR band is only lightly affected by extinction. In Fig. [9] (bot- 
tom panel) we plot the monochromatic 6 fim luminosity (Le /^m) 
of the AGN, derived from the SED fits, as a function of the X-ray 
luminosity (2 - 10 keV, rest-frame; not corrected for absorption). 
Where the AGN component was not significantly detected in the 
IR band, an upper limit was calculated for to be 30% of the 
total (AGN + SFG) 6 fim luminosity estimated from the best-fit 
SEDs (Table [T[ see Appendix A for justification). The shaded 
area in the figure represents the typical X-ray/IR luminosity cor- 
relation found for local, unobscured AGN by |Lutz eT al. (2004), 
while the dashed line represents the Lf, ;jm-Lx relation predicted 
for sources obscured by large column densities (A^h ~ 10^"* cm"^; 
[Alexander et al.|2008 1. When sources are heavily obscured (e.g., 
Compton thick), their X-ray luminosity can be strongly sup- 
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Fig. 9: Fraction of IR detected AGN as a function of the X-ray luminosity (top panel); the sources have been divided in five luminosity 
bins, between Lx = 10'*" - lO'*^ erg s"', although no AGN is expected to be found below Lx = lO'*' erg s"'. The bottom panel shows 
the X-ray luminosity at 2 - 10 keV (rest-frame; not corrected for absorption) vs. the monochromatic 6 /urn AGN luminosity (vLv) for 
the radio-excess sources. Some local Compton-thick AGN are also plotted as open squares (thin marks: observed Lx; thick marks: 
absorption corrected Lx). The shaded region represents the typical L(, - Lx relation for local AGN ( Lutz et al.|2004 i, while the 
dashed line is the relation expected for heavily obscured AGN (A^h ~ 10^^ cm"^) and the dash-dotted line is the typical relation for 
starburst galaxies (taken from Alexander et al. 2008[l. 



pressed, as opposed to the IR luminosity, and they should lie on 
the left hand side of the local AGN relation (or even on the left 
of the dashed line). Indeed, for comparison, we also show in the 
figure some local Compton-thick AGN, with observed (thin open 
squares) and intrinsic (i.e., corrected for absorption; thick open 
squares) X-ray luminosities, taken from literature (NED''). Since 
for these local CT AGN we did not perform SED decomposition, 
we plotted the total 6 yum luminositjj^for these sources (obtained 
from Spitzer or ISO data, taken from literature). 

About 59% of our radio-excess AGN (19 out of 32 source^ 
plotted in Fig. |9) bottom) follow the local Le ^/m-Lx relation from 
[Lutz et ar (2004), indicating good agreement between the X-ray 
luminosity and the MIR luminosity estimated from our SEDs; the 
remaining X-ray detected AGN (8 sources) lie closer to the re- 
gion expected for heavily obscured AGN (A^h ~ lO^'* cm"^; see 



* For sources like Mkn 3 and Mkn 23 1 , this is a good estimate of the 
AGN 6 fim luminosity since their MIR spectra are dominated by the 
AGN (e.g. Gouldin g eraLl[20T2l [Redman et al.|2005| l; however, for 
other sources, e.g. NGC 4945, where the MiR emission is dominated 
by star formation (e.g. [Goulding et al.|[Wl2, Peeters et al.||2004l > the 
plotted 6 jim luminosity is likely to be an overestimate of the intrinsic 
AGN emission at 6 yum. 

' The 32 sources plotted in Fig. |9]include all of the radio-excess AGN 
with a measured AGN luminosity at 6 /jm and/or in the X-ray band; 
namely, they are 23 IR AGN, 7 X-ray detected IR SFGs and 2 X-ray 
detected passive systems (see Table|2j. 



Alexander et al. 2008[ l, together with all the X-ray undetected 
sources (5 sources), for which we have X-ray luminosity upper 
limits. The MIR luminosity of these objects indicates that the 
intrinsic power of the AGN is higher than that detected in the 
X-rays, suggesting they may be heavily obscured, possibly by 
Compton-thick material surrounding the black hole (A^h > lO^'* 
cm.-^). Also for the X-ray undetected IR AGN (-22%, 5/23 
sources; see Table|2]i the IR luminosity measured from the SEDs 
(vLfj^n, > 10'*^ erg s"') implies that these sources should be 
bright enough in the X-rays to be detected by the deep Chandra 
data available in GOODS-N. Since our capability to detect the 
AGN emission in the IR decreases with decreasing AGN power 
(Fig. |9] top), this suggests that the five IR AGN that are X-ray 
undetected are likely to be intrinsically luminous sources, and 
therefore they are candidate Compton-thick AGN. 

On the basis of these analyses, amongst the IR detected AGN 
we can therefore estimate the fraction of candidate Compton- 
thick AGN to be up to ~43% (10 sources that lie close to the 
A^H ~ lO^'* cm"- line in Fig. [9] bottom), considering both 
X-ray detected and undetected AGN. However, to confirm the 
classification of these AGN as Compton-thick and their frac- 
tion, other diagnostics are necessary (e.g. Georgantopoulos et al. 
|201 1[ [Goulding et al.|[2bll| l. The most efficient way to unam- 
biguously classify these sources as Compton-thick AGN, would 
be the detection of a strong iron line (Fe Ka, E - 6.4 keV) and a 



Article number, page 13 of 



28 



reflection dominated spectrum from detailed X-ray spectral anal- 
ysis. Unfortunately, since our sources are at high redshifts, they 
are typically weakly detected (or undetected) in the deep Chan- 
dra data, preventing the good quality spectra necessary for this 
kind of analysis. 

5.2. Infrared star-forming galaxies 

Together with the IR AGN, the other main group amongst our 
radio-excess sources {^45%) have IR SEDs consistent with star 
formation emission, with no significant evidence at IR wave- 
lengths for the presence of an AGN. However, the excess ra- 
dio emission still suggests that these sources are not simple star- 
forming galaxies. 

|Miller & Owen ( [200 1') found a significant number of star- 
forming galaxies not hosting AGN (~ 20%) with enhanced radio 
emission compared to the typical FIR-radio relation. However 
this phenomenon has been observed in the centre of galaxy clus- 
ters and the radio excess for these galaxies has been estimated 
to be not larger than a factor of three from the field galaxy FIR- 
radio relation (corresponding to <7 ^ 1.73; |Gavazzi & Jaffe 1986). 
[Bressan et al.| (|2002 ) interpret this efifect as a post starburst phase, 
when the star formation has stopped, while the radio emission 
produced by supernova remnants (SNRs) has still not dimmed, 
as it fades more slowly. This phenomenon can not explain the 
origin of the radio excess in all of our IR SFGs, since the selec- 
tion criteria adopted in our analysis (^ < 1 .68), which correspond 
to more than a factor three enhancement of the radio emission, 
should already exclude the majority of these radio-excess post- 
starburst galaxies; moreover, it is unlikely that about half of our 
radio-excess sources reside in the centre of galaxy clusters. In 
fact, the average q calculated for our IR SFG ((^)ir sfg = 1.33; 
see Sect. |4.3[ ) suggests that the radio emission in these sources is 
much higher than that observed in post-starburst galaxies. More- 
over, we verified that the SFR for the IR SFGs is consistent with 



that of star-forming/starburst galaxies (Sect. 4.3 i, meaning that 
the star formation is still ongoing in these systems. Therefore, 
the excess radio emission in these IR SFG can only be due to the 
presence on an AGN. 

A number of the IR SFGs (7/23) are detected in the X-ray 
band, although they are typically less luminous than the X-ray 
detected IR AGN (Lx ^ 10'*^ erg s"'). For these sources we es- 
timated the IR luminosity expected for the AGN at 6 fim from 
the X-ray luminosity in the 2-10 keV band (rest-frame) using 
the fLutz et al.| ( |2004j ) relation for local unobscured AGN. In four 
cases, the estimated Lg of the AGN is consistent with the total 
6 fim luminosity calculated from the SEDs (see Fig. |9]l, indi- 
cating that the MIR emission of these sources is unlikely to be 
produced by star formation only, despite our SED fits suggest 
that this is the case. The AGN in these sources might be ob- 
scured at MIR wavelengths (e.g. Deo et aL]|2009 |Goulding &| 
I Alexander||2009 [ |Goulding et al. |20I2| i and therefore the emis- 
sion appears dominated by star formation. These IR SFGs are 
the brightest sources amongst this grou J"] (Lx ~ lO'*-' erg s"') 
and indeed their SEDs could be consistent with a combination 
of SFG H- AGN emission; however, given our conservative crite- 
ria to identify the best-fitting SED model, the AGN component 
in these cases did not pass the /-test probability threshold (see 
Appendix A). At the luminosities of these sources, the SED fit- 
ting procedure struggles to disentangle the two components, es- 
pecially if the AGN contribution is not significantly larger than 
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We stress that the X-ray luminosities are not corrected for absorption 
(Sect. |2.4| l, and therefore if large amounts of dust and gas are present in 
these sources, the intrinsic X-ray luminosity is likely to be higher. 
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Fig. 10: Specific star formation rate (sSFR) vs. redshift for the 
entire VLA/24 /im sample (grey dots) and for the radio-excess 
sources, divided into their three IR categories: IR AGN (red 
squares), IR SFG (blue circles) and passive systems (asterisks). 
The four sources plotted at the bottom (with a symbolic value 
of sSFR= 0.001 Gy-i) are the four IR AGN for which the SFR 
could not be estimated (see Sect. |4.3[ ). The tracks in the plot 
show the sSFR evolution with redshift for main sequence star- 
forming galaxies (MS; solid line) and for star-forming galaxies 
(SFG; dotted line) from |Elbaz et"aL] ( |201 Ij ). The radio-excess 
sources have typically lower star formation than the star-forming 
galaxy population detected at 24 yum. 



that from star formation (Fig. [9] top, and Appendix A). For the 
other three sources, both the X-ray and the predicted IR lumi- 
nosities are below Lx 5 10"*^ erg s"' (see Table [l]), which are 
consistent both with a faint AGN or with star-forming galaxy 
emission. The remaining »70% of the IR SFGs (16/23) are not 
detected in the X-ray band. 

The most likely interpretation for these IR SFGs amongst our 
radio-excess sources is that they are low-moderate luminosity 
AGN, typically with Lx < lO^^* erg s ' (see Fig. |9]l, or possi- 
bly luminous AGN heavily obscured at MIR wavelengths, and 
hosted in dusty star-forming galaxies, whose IR emission is thus 
overwhelmed by star formation and eludes our detection capa- 
bility (Sect. |5.1| l. If a large amount of dust is present in these 
galaxies, the nuclear emission from the AGN could be heavily 
obscured and remain undetected even in the X-rays. Supporting 
this interpretation is the fact that one of our X-ray undetected IR 
SFGs (#29, in Table [TJ was found to have a compact radio core 
in the VLBI observations by Chi et al. ( 2009[ ), clearly revealing 
the presence of an AGN in this star-forming galaxy (see Sect. 
|4.2| l; this source also shows extremely red optical/NIR colours, 
indicating significant dust reddening ( |Chi et al.|2009[ l. Although 
the inability to unambiguously identify the presence of an AGN 
in the IR and/or X-ray bands for the majority of these IR SFGs 
represents an evident limitation in finding AGN, the radio excess 
provides a clear signal, often the only signal, of nuclear activity 
in these sources, bringing us closer to completing the census of 
the AGN population. 
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5.3. "Passive" systems 

A small fraction (a!lO%, 5/51 objects) of the radio-excess sources 
seem to have weak contributions from AGN activity and star for- 
mation, as their IR SEDs are in fact consistent with those of rel- 
atively quiescent elliptical galaxies (Sect. 



passive from our SED fitting analysis have, in fact, sSFR< 10% 



4.3 I 



These sources 

have lower redshifts (z < 1, except for one with z ~ 1.5) than 
the average of the radio-excess sample and some of them are 
amongst the most powerful radio emitters (Table 1). This sug- 
gests that they might be radio galaxies, such as Fanaroff-Riley 
type I radio galaxies (FR Is), which are powered by low accretion 
rate AGNs ( |Chiaberge etal.|1999|piisellini & Celotti|2001| i. At 
low accretion rates (L/LEdd < 0.01) AGNs are predicted to have 
radiatively inefficient accretion flows (RIAFs; 'Narayan & Mc- 
Clintock|2008 Narayan et al.|1995 1, which are typically associ- 
ated with st rong radio outflows ( Maccarone et al. 2003 ; Narayan 
et al. II 19951). Indeed, one of the passive sources in our sample 



(#47 in Table [TJ is a well known wide-angle-tail (WAT) radio 
source at z = 1.265 and shows strong extended radio jets (e.g., 
Bauer et al.|2002[ ). These galaxies may have already passed the 
active phase when the SMBH is rapidly growing and the host 
galaxy experiences strong star formation and they are left with a 
low accreting SMBH in a passive red galaxy (e.g. Hickox et al. 
2009] [Hopkins et al.|2006a| ). 



5.4. Radio excess or dimmed star formation? 

In this subsection we aim to investigate the properties of the 
galaxies hosting the radio-excess AGN, comparing them to the 
properties of X-ray selected AGN hosts. In Section 4.3 we have 
shown that the radio-excess sources have on average star forma- 
tion rates typical of normal star-forming/starburst galaxies. How- 
ever, to have better insights into the host galaxy properties one 
should explore the specific star formation rate (sSFR), i.e. the 
SFR per unit stellar mass. In fact, recent studies have shown that 
a correlation exists between the SFR and the galaxy stellar mass 
(e.g. [Daddi et al.|[2007l pbaz et al.|[2007l |Noeske et al.|[2007l 
Pannella et al.|2009a i, meaning that more massive galaxies have 
typically higher SFRs. This correlation also evolves with red- 
shift, producing the so called "main sequence" of star-forming 
galaxies (MS; e.g., |Elbaz et al.|201 1 ; Noes ke et al.|2007[ ). There- 
fore, removing the mass dependence of the SFRs gives a better 
estimate on the relative star formation activity of galaxies. 

We calculated the sSFR by simply dividing the SFRs (Sect. 

23] and Fig. [T) 



4.3 I by the galaxy stellar masses (see Sect. 



Figure 10 we show the sSFR as a function of redshift for the 
entire VLA/24 fim sample (grey dots) and for the radio-excess 
sample, distinguishing between the different types of sources ac- 
cording to their IR SEDs (i.e., IR SFG, IR AGN and passive). 
In the plot we also show the sSFR-z evolution track for MS star- 
forming galaxies and for starburst (SB) galaxies defined by Elbaz 



et al. (2011 1. Our radio-excess sources typically lie below the 



sSFR tracks for star-forming galaxies, indicating that they have, 
on average, lower star formation contribution when compared to 
radio-normal sources of the same stellar mass. However, if we 
consider the average sSFR for the IR SFGs, IR AGN and passive 
systems amongst our radio-excess sample, we find that for the IR 
SFGs <sSFR)iR SFG = 1-15 + 0.26 Gy"', consistent with the MS 
star-forming galaxies (assuming an average redshift of z » 1.5), 
while for the IR AGN, the average sSFR is significantly lower: 
<sSFR)iR AGN = 0.38 ± 0.14 Gy"'[3 The sources identified as 



' ' The four IR AGN where the SFR could not be calculated (see Sect. 
[4.3^ were not included in the average sSFR calculations, and are plotted 
in Fig. fiolwith a symbolic value of sSFR= 0.001 Gy"'. 



of the average MS galaxies at a given redshift (Fig. 10 1 and are 
more massive ((M,) a; 11.5 M©; see Table [TJ than the rest of the 
sample «M.) a; 11.0 Mo and <M.) a 10.7 for IR AGN and 
IR SFGs, respectively)!^ 

Previous studies have found that galaxies hosting X-ray se- 
lected AGN typically have star formation rates in agreement with 
the MS star-forming galaxies (e.g. [Alonso-Herrero et al. 2008 



MuUaney et al. 



20T2||Santini et al. "2012; Si lverman et'aLi2009 



Xue et al.|2010 and do not show different properties from galax- 
ies not hosting AGN. This suggests that the radio-excess tech- 
nique tends to select a different AGN population from those se- 
lected in X-rays, with the host galaxies of radio-excess AGN 
having on average less ongoing star formation when compared 
to the hosts of X-ray selected AGN. We note that thanks to our 
SED decomposition, we were able to calculate the SFRs (and sS- 
FRs) more accurately than in many of the previous works men- 
tioned above by removing any AGN contribution to the IR lu- 
minosity (see Sect. |4.3| l. This gives us smaller sSFR values 
than those resulting assuming a pure SFG dominated IR lumi- 
nosity. However, we found that the discrepancies between the 
sSFR of radio-excess AGN and X-ray AGN hosts are real and 
are not due to the technique adopted to calculate the L/^ (and 
therefore sSFR). In fact, if we estimate the sSFRs for all of the 
X-ray detected sources in the VLA/24 fim sample (see Sect. |2.4| l, 
we obtain good agreement with the MS sSFR (Del Moro et al., 
in prep.). We therefore conclude that the host galaxies of radio- 
excess AGN are growing at a slower rate than the typical X-ray 
selected AGN hosts. Since many of our radio-excess sources are 
not detected in X-rays, it seems that the excess radio emission 
can reveal a complementary AGN population to that selected in 
the X-ray band. 

Since in the local Universe more radio-loud AGN are pref- 



erentially found in galaxy groups and clusters (e.g. Best||2004 



jKauffmann et al.|[2008| ), as opposed to typical radio-quiet AGN 
that tend to reside in less dense regions (e.g. KaufFmann et al.j 



2004 1, we want to verify whether the radio-excess AGN also 



reside in different environments than the typical X-ray selected 
AGNs. We therefore cross-matched the positions and redshifts of 
our radio-excess sources with those of extended X-ray sources, 
coiTesponding to known clusters, in the GOODS-N field (e.g. 
jBauer et al.|2002| l. As search regions we used the ellipses defin 



ing the extent of the X-ray emission, reported in Table 1 of Bauer 
jet al.| ( 2002| . We found that only two out of the 51 radio-excess 



In AGN in our sample (specifically sources #6 and #17 from Table 



lie within the X-ray extended emission regions from Bauer] 
et al.| ( ]2002] l; only one of these two sources (#17; i.e. only ^2% 



of the radio-excess sample) has a redshift matching that of the 
galaxy cluster (z = 1.01; Bauer et al. 2002] l and therefore is likely 
to be associated with the cluster. The large majority of our radio- 
excess sources lie in less dense environments, similarly to the 
X-ray selected AGN population. 

The lower sSFRs observed for our sources compared to MS 
galaxies and X-ray AGN hosts suggest that the radio-excess char- 
acteristic of our sources is not only due to a stronger radio core 
emission from the AGN than in "radio-normal" sources, but also 
to a dimming of the emission from star formation, which en- 
hances, by contrast, the radio excess. The decreasing average 
sSFR for the different types of radio-excess sources identified 
through our SED analysis (i.e. IR SFG, IR AGN, passive; Sect. 



We note that the range of stellar mass values for the radio-excess 
sources is consistent with that of the entire VLA/24 yum sample (see 
Sect. |2.6| l; however, on average, the radio-excess sources have larger 
stellar masses (median M, a 11.1 Mq) than the "radio-normal" sources. 
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|4.3| l indicates that we might be looking at different stages of 
BH-galaxy evolution: from the IR SFG that are the more ac- 
tive and possibly the most dust-obscured sources (see Sect. |5.2| i, 
where the presence of the AGN is testified only by the radio- 
excess, to the IR AGN, where the dust obscuration is progres- 
sively lower (although still very high in some cases; see Sect. 
|5.1| l, and the AGN is detected at radio, IR and often X-ray bands, 
while the star formation is quenching. At the final stage there are 
the passive systems, where the star formation has stopped and 
the SMBH is slowly accreting in a massive, passive galaxy (e.g. 
Hickox et al.|2009| [Hopkins et al.|2006a"l l. 



6. Summary and conclusions 

Using deep IR Spitzer and Herschel data, as part of the GOODS- 
Herschel program, and radio VLA data in the GOODS-N field, 
we performed a detailed SED analysis for 458 sources with 
spectroscopic or photometric redshift identification (z < 3.0). 
From the best-fit SEDs we calculated the FIR flux (/fir; A - 
42.5 - 122.5 jum) and the FIR-radio relation {q) to identify a 
sample of sources with excess radio 1 .4 GHz emission over that 
expected for star- forming galaxies {q « 2.2; 'Helou et al. 1985 ). 
We obtained a sample of 51 (^11% of the initial sample) radio- 
excess AGN (<7 < 1 .68) and investigated their radio, IR and X-ray 
properties. The main results of our analysis can be summarised 
as follows: 

- We found that the fraction of radio-excess sources increases 
with X-ray luminosity (Sect. |4.2| i, suggesting that the more 
luminous AGN are also more powerful in the radio band, al- 
though not necessarily radio-loud. This suggests that there 
is a wide distribution of radio power amongst the AGN pop- 
ulation and that the radio-excess sources (i.e. sources with 
intermediate radio power) constitute a significant part of the 
total AGN population. 

- The radio-excess sample seems to be composed of a hetero- 
geneous mix of sources: i) IR AGN (^A5%, 23/51), which 
are the "classical" highly accreting AGN hosted in star- 
forming galaxies; ii) IR SFG (~45%, 23/51), which are likely 
to host moderate luminosity (Lx < 10"*^ erg s"') or dust- 
obscured luminous AGN, whose emission eludes our detec- 
tion capabilities in the IR and often the X-ray bands, but they 
are identifiable through their radio excess; iii) low accretion 
rate AGNs (e.g. RIAFs, FR Is) hosted in passive, non-star 
forming galaxies (^10%, 5/51), where the excess radio emis- 
sion is likely to be due to the presence of strong radio jets and 
lobes, in addition to a radio core (Sect. [5]l. 

- Only 27 of the radio-excess AGN (~53%) are detected in the 
X-ray band; the large number of X-ray undetected sources 
suggests that many radio-excess sources might be heavily ob- 
scured (possibly Compton-thick) AGN (Sect. \5.\\ . Indeed, 
amongst the 24 X-ray undetected sources, five sources (i.e. 
~20%) show a strong AGN component in the IR band. The 
majority of the X-ray undetected sources are IR SFGs, which 
are likely to be low-moderate luminosity AGN, or luminous 
dust-obscured AGN, as demonstrated by the detection of a 
compact radio core in deep VLBI observations for one of 
these X-ray undetected IR SFGs (Sect. |5.2| l. Assuming that 
the local correlation between the AGN Le^jm and the X-ray 
(2-10 keV) luminosity holds out to high redshift, we esti- 
mate the fraction of Compton-thick candidates amongst our 
radio-excess AGN to be ~20% (~43% amongst the IR de- 
tected AGN 10/23); however the unambiguous classification 
of these sources requires other diagnostics, which are not fea- 
sible with the current data (Sect. |5.1| l. 



- The specific star formation rates (sSFR) estimated for the 
radio-excess AGN are on average lower than those observed 
for main sequence star-forming galaxies, or X-ray AGN host 
galaxies, indicating that the radio-excess technique tends to 
select a different AGN population than the X-rays (Sect. 5.4 1. 
This also suggests that the excess radio emission measured 
for our sources is not only due to strong radio emission from 
the AGN, but to a combination of two effects: i) radio core 
emission from the AGN and ii) dimmed star formation emis- 
sion. Moreover, the progressively lower sSFR observed for 
IR SFGs, IR AGN and passive systems, suggests that through 
our radio-excess sources we might look at different stages of 
BH-galaxy co-evolution. 
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Appendix A: Testing tlie SED fitting technique 

As described in Section [TT] our fitting approach consists of us- 
ing 5 SFG templates, derived from a sample of local star-forming 
galaxies and an empirically derived AGN template ( |Mullaney| 
et al.||2011 ) modified by an extinction law (with Ay ~ 0-30 
mag). We applied these templates to decompose the total SEDs 
of high redshift sources (z < 3) into AGN and galaxy compo- 
nents. To determine the best-fitting solutions, we firstly fit to the 
Spitzer and Herschel photometric data a simple model including 
only the SFG template, using minimization, yielding 5 SED 
solutions (one for each SFG template). We then fit the data with 
a more complex SFG + AGN (plus extinction on the AGN com- 
ponent only) model, using an /-test to determine the best-fitting 
solutions. The criteria we used to choose whether the AGN com- 
ponent significantly improves the fit are: i) an /-test probability 
>90% confidence when adding the AGN template to the fit; ii) 
this first condition must be met in the majority of the SED fitting 
solutions (i.e. at least 3 out of 5). 

This SED fitting approach (Sect. |3.1| l and the criteria cho- 
sen to identify the best-fit models were selected based on the 
results of careful tests performed on a sample of 30 sources in 
the VLA/24 detected sample (Sect. [2]) for which Spitzer 
IRS low-resolution spectroscopy is available (covering the wave- 
length range /I a; 3 - 20 jum, rest-frame; Kirkpatrick et al.|2012[ 
Murphy etalT2 009; Pop e et al.|2008|). These sources are all de- 
tected in the X-ray band ( {Alexander et al.|2003| l and the majority 
of them are also detected in at least one of the Herschel bands 
at 100, 160 and/or 250 ^im (29/30 sources; Table [Al]). Our tests 
consist of two main steps: 

1. constraining the AGN and star formation components using 
the IRS spectra (plus Herschel data), which provide more de- 
tailed information (e.g. PAH features or featureless power- 
law continuum) than photometry alone on the source emis- 
sion at MIR wavelengths, where the AGN contribution af- 
fects the SED the most (/I < 30 - 40 yum); 

2. performing the SED fits using only Spitzer and Herschel pho- 
tometric points (not including the IRS spectra), testing var- 
ious criteria to select, amongst the different solutions, the 
best-fitting models that more closely matched the IRS spec- 
tral fitting results. 

In analysing the IRS spectra, we also included the Herschel 100, 
160 and 250 //m flux densities to help constrain the SEDs in 
the FIR band, as the spectra alone cover too small a wavelength 
range to allow a reliable extrapolation of the SEDs in the broad 
IR band. In fact, fitting just the IRS spectra (not including the 
Herschel data) tend to underestimate the amount of star forma- 
tion contribution to the total emission, favoring higher contribu- 
tion from the AGN component, instead. 

Since the Spitzer IRS spectra provide a large amount of data 
for fitting our SED templates, we fitted the spectra (and Herschel 
data) directly using a model including both SFG and AGN (plus 
extinction) templates described in Sect. |3.1| to measure galaxy 
and AGN contributions to the total IR emission. We note that be- 
cause the IRS data have much smaller uncertainties compared to 
the Herschel data, we increased the errors on the IRS spectra by 
a constant factor to be of the same order of magnitude of the Her- 
schel data error ^(when S/N> 3); this is to avoid the SED fitting 
results being dominated by the Spitzer IRS data (which would 
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Fig. A.l: Comparison of the AGN 6 /jm luminosity obtained 
from the SED fitting of the photometric Spitzer and Herschel 
data points (vLe /jm, SED) versus the L(, estimated from the 
IRS spectral (-1- Herschel data points) fitting (vLe fim, IRS). For 
the IR SFGs, where the AGN component was not significantly 
detected, we plotted an upper limit of the AGN Lg ^n, (i.e., 30% 
of the total 6 fj.m luminosity; see Appendix A). Sources where 
the classification from the two analyses agree are shown as filled 
symbols: IR AGN (filled squares) and IR SFG (filled circles). 
The AGN identified only from the IRS spectral fit are shown as 
open squares. 



have much more "weight" on the resulting than the Herschel 
data), and also to account for the intrinsic scatter between our 5 
discrete SFG templates. We adopted minimisation to evaluate 
the best fit of the 5 SFG templates. We note that due to the high 
signal-to-noise ratio of the IRS spectra (despite the enhancement 
applied to the errors), the resulting values for all the fits are 
very high (x^/d.o.f. » 1) and therefore cannot be used as an 
abs olute measure of the go odness of the fit in the usual way (see 
also 



MuUaney et al.|20Tr i. However, since the;^'^ calculation is 



'^^ For the Herschel undetected source, i.e. where S/N< 3 in the Her- 
schel bands, the errors on the 100, 160, 250 /im flux density measure- 
ments are still much larger than the Spitzer IRS spectra "increased" un- 
certainties. 



affected by this issue in the same way for each solution, a com- 
parison between x^ values can still be used to identify the best-fit 
amongst the five different solutions: the spectral fit yielding the 
minimum x^ was chosen as the best-fitting solution. It is im- 
portant to note however, that in some cases the difference in x^ 
between different solutions is small, and therefore, the chosen 
best-fit model is not the only acceptable solution to reproduce 
the data. To overcome this issue, instead of choosing one solu- 
tion as the best-fit, we used a weighted mean between the five 
different fitting solutions to retrieved the 6 jim luminosity of the 
AGN (from the fitted AGN component), the FIR fluxes and the 
q values for the sources. The;\f^ values were used to calculate 
appropriate weights in order to give the maximum weight to the 
solution with the minimum (i.e. u;,- = X^mi„lx])- ^ dominant 
AGN component (i.e. Lg /jm,AGN/i'6 //m.tot > 0.5) was measured in 
16 cases. From these IRS spectra (-^-Herschel) fits, we also deter- 
mined that when the AGN is found to contribute less than ~3Q% 
to the total (AGN -h SFG) emission at 6 yum, the total SED is 
not significantly aff'ected by the presence of the AGN component 
compared to that of a pure SFG. Therefore, when the AGN is not 
significantly detected with our SED fitting approach (i.e. when 
using photometric Spitzer and Herschel data only), we chose the 
upper limit for the AGN luminosity at 6 pm to be 30% of the 
total 6 pm luminosity (see Sect. 5.1 1. 
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After using the IRS spectra to obtain good constraints on the 
AGN and star formation components, we then performed sev- 
eral sets of SED fits to these sources using Spitzer 8, 16, 24 jum 
and Herschel 100, 160, 250 yum flux densities only (not including 
the IRS spectra), testing different criteria to select the best-fitting 
solutions, in particular to establish the significance of the AGN 
component in the fits. These tests were performed to find the best 
"recipe" to recover, with these sparse photometric data points 
(which is the approach adopted in this paper; see Sect. |3.1| l, the 
results obtained from the IRS spectral (plus Herschel data) fits 
(Fig. 



To test the results for our chosen SED fitting approach, we 
compared the 6 /vm luminosity of the AGN obtained from the IRS 
spectral fitting (+ Herschel data) to that obtained from the SED 
fitting performed using Spitzer and Herschel photometry (Figure 
A.l I. In Fig. A.l we distinguish between sources with a domi- 
nant AGN component at 6 //m (IR AGN: Le nm, agn/L(, fim, tot > 
0.5) from those where the AGN component is not significantly 



detected (IR SFG: Le/im, AGt^/L(, 



fim, tot 



< 0.5). In a;87% of the 



A.2 and Table A.l 



As for the IRS spectra, the Spitzer 
photometric data have much smaller uncertainties than the Her- 
schel data; we therefore increased the errors on the 8, 16 and 24 
fj.m flux densities (by constant factors for all of the sources) dur- 
ing the SED fitting process in order for each data point to have 
a similar "weight" in the fit and to allow for the intrinsic scat- 
ter on the templates. We note, however, that in these SED fits, 
as for the IRS + Herschel data fits, the reduced alone (i.e. 
X^/d.o.f. » I) is not a reliable indicator of the goodness of fit 
since the is very sensitive to the errors on the data points used 
in the fits, especially when the number of data points is limited, 
as in our case. Nevertheless, as we pointed out above, the relative 
can still be used to determine the best-fit model between two 
sets of solutions (SFG only, or SFG + AGN; Figure [3]), so we 
performed an /-test, which compares the;^'^ values and d.o.f. of 
the two models (see Sect. 



cases the IR classification of the sources obtained from the IRS 
spectra and the SED fits are consistent with each other (12 IR 
AGN and 14 IR SFGs) and the 6 ijm luminosities measured from 
the two types of analysis are in good agreement (Table [AT| . In 
the remaining cases (a; 13%), the two different analyses yielded 
different results: the IRS spectral analysis identified an AGN 
component in four further sources, which are missed by our SED 
fitting analysis; this is because the criteria we adopted in defining 
the AGN component as significant in our SED fits are very con- 
servative and thus, while we can easily identify the AGN where 
Le fim, AGt^ / L(, fjm, tot » 0.5, wc are likely to miss some cases 
where the AGN component is not strongly dominant in the MIR 
band. With this in mind, we can consider the detection of the 
AGN emission component through our SED fitting approach as 
robust. 

As a final verification of the reliability of our SED fitting ap- 



3.1 



proach, we plot in Figure A.2 the best-fit SEDs (using Spitzer and 



to measure the improvement of 
the fit obtained by including the AGN component to the model. 

We tested the reliability of the best-fitting solutions by vary- 
ing the /-test probability threshold between the two adopted 
models (SFG and SFG -i- AGN) and also varying the number 
of solutions that needed to pass this threshold for the model to be 
accepted as the best-fit. A higher confidence level (e.g. 95% or 
99%) allowed us to recover the AGN component detected in the 
IRS spectra only in a small number of cases, while missing some 
very obvious, strong IR AGN, where the AGN emission found 
from the IRS contributes 100% to the total 6 fim luminosity. We 
therefore lowered the confidence level threshold from the /-test 
to be at least 90%. The /-test was performed between each pair 
of SED fitting solutions (i.e. for each SFG template) in order 
to verify whether the need of the AGN component was depen- 
dent on the SFG template used in the fits. Indeed, we recognised 
that the fits performed with one particular SFG template ("SB5" 
in |Mullaney et al.|[20II[ ) typically required an AGN component 
with higher confidence than the others. To avoid biases on the 
results due to the choice of SFG template and thus to remove 
the dependence of the fitting solutions on any specific templates, 
we decided to accept the AGN component as significant only if 
the /-test criterion was met in at least 3 of the 5 SED fitting 
solutions, otherwise we conservatively assumed the simple SFG 
model as the best-fit. 

As previously seen from the Spitzer IRS spectral fits, in 
some cases different fitting solutions yielded small differences 
inx^ values meaning that a unique solution could not be defined. 
Therefore, once we established the best-fit model (SFG or SFG 
+ AGN), we obtained the 6 pm luminosity of the AGN, the to- 
tal FIR fluxes and the q values for our sources using a weighted 
average of the values calculated from all the best-fit model solu- 

tion^(as for the fits of the Spitzer IRS spectra + Herschel data); we only plotted the SED with the lowest x^ value amongst the 



Herschel photometric points only) for the 30 sources analysed 
here in comparison to the Spitzer IRS spectra; we showed only 
one of the best-fitting solutions for each source. The SEDs were 
fitted using only Spitzer 8, 16, 24 fim, and Herschel 100, 160, 
250 pm photometric points (black circles); the Spitzer-IRAC data 
points at 3.6, 4.5, and 5.8 pm (red triangles) and the VLA radio 
data points (black star) are over-plotted to the SEDs (not used 
in the fits) to show the agreement with the resulting SEDs. The 
Spitzer IRS spectra (cyan line) are also shown, but they are not 
included in the SED fits. The panels on the top right of each plot 
is a zoom on the IRS spectra (2.5 - 25 pm rest-frame) to better 
show the comparison between the best-fit SEDs and the spectra. 
In general, there is very good agreement with the resulting SED 
and the IRS spectra, even though we stress that here the spectra 
were not used to constrain the SEDs. This result is a confirma- 
tion of the validity of our templates and SED fitting approach. 
The only few cases (2/30; i.e. CXOJ123555. 13+620901.7 and 
CXOJ123726.51-H622026.8) where the SEDs diff^er from the 
MIR IRS spectra are those where deep silicate features are 
present; in these cases the photometric data points do not provide 
enough information to predict the amount of extinction needed to 
reproduce these strong features. However, since the purpose of 
our analysis is not to measure the strength of the spectral fea- 
tures, such as silicate absorption/emission features or PAH emis- 
sion lines, but to estimate the AGN and SFG contribution to the 
overall IR SEDs, we can ignore the discrepancies between SEDs 
and spectra for these sources. 

Appendix B: SEDs of the radio-excess sources 

In this appendix the best-fit SED plots for the entire radio-excess 



sample (51 sources) are reported (Figure B.l i; for each source 



the errors on the estimated averages were calculated as weighted 
average variances. 



best-fitting model solutions (see Sects. |3.1|and Appendix A) 



For the SFG + AGN model we only included in the weighted average 
calculation the solutions where the AGN component was significant in 
the fit (>90% confidence level from the /-test). 
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A. Del Moro et al.: Radio-excess sources in GOODS-Herschel. 
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Fig. A. 2: Best-fitting solutions to the spectral energy distributions (SEDs) of the X-ray detected sources in GOODS-N with Spitzer- 
IRS spectroscopy (Kirkpatrick et al. 2012; Murphy et al. 2009 ; Pope et al.|2008 1. The dotted lines represent the AGN component 
and the dashed lines indicate the SFG component; the total SEDs are represented as solid lines. The Spitzer 8, 16, 24 yum and the 
Herschel 100, 160, 250 pm data points (black circles) have been used to constrain the SEDs, while the Spitzer-IRAC 3.6, 4.5, 5.8 
yum (red stars) and the VLA 1.4 GHz (black star) radio data points are over plotted on the SEDs but they are not included in the fit. 
The Spitzer-l^S spectrum is also shown (cyan line), but it is NOT used to constrain the SEDs. On the top right-hand side of each 
plot, a zoom on the Spitzer-]RS spectrum is shown (2.5-25 pva); our best-fit SEDs are typically in very good agreement with the 
Spitzer-IKS spectra. 
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Fig. A. 2: Continued. 
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Fig. A. 2: Continued. 
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Fig. B.l: Best-fit SEDs for the 51 sources in the radio-excess sample. The total SEDs are shown as black solid lines, the AGN 
templates are shown as dotted lines and the SPG templates as dashed lines. Filled circles represent the Spitzer 8, 16, 24 jim and 
the Herschel 100, 160, 250 /im flux densities, which are used to constrain the SEDs. Open symbols indicate the data that were not 
included in the SED fitting process: red triangles are Spitzer-IKAC 3.6, 4.5, 5.8 //m flux densities, black open circles are SPIRE 350 
and 500 yum, and black squares are VLA 1 .4 GHz flux densities; the blue stars represent the 6 /im luminosity of the AGN predicted 



from the X-ray luminosity (2-10 keV, rest-frame) using the Lutz et al. (2004) relation for local unobscured AGN; we note that these 
points do not always match the IR AGN component because the X-ray luminosity tends to underestimate the intrinsic AGN power if 
the AGN emission is heavily absorbed. For the passive sources we plotted the SFG template upper limit (grey line) and an elliptical 
PL'rtM^Afifii&^f|'-^a^&5§ ^IfP^'^ 4-3 1 to show that it could well represent the data, although we stress that this latter template 

is not includeci in our SETB-fitting analysis. The rise of the IRAC data point at short wavelengths observed for many sources is due to 
the emission from the galaxy old stellar population. 
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Fig. B.l: Continued. 
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Fig. B.l: Continued. 
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